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ABSTRACT 

This  thesis  describes  the  design  and  construction  of  a  wide  band 
(60-3060  Hz),  three  phase  (including  zero  sequence),  non-transposed, 
lumped  element,  power  transmission  line  model.  The  model  consists  of 
from  15  to  65  "three  phase  L  sections"  -  depending  upon  the  length  of 
the  actual  line  being  modeled  -  and  includes  trimming  features  to  model 
both  horizontal  and  vertical  configurations  of  lines  in  the  345  kV  to 
1500  kV  range. 

To  predict  the  performance  of  the  model   and  determine  the  number 
of  sections  required  to  model   a  given  line,  the  poles  of  the  driving 
point  impedance  of  various  numbers  of  L-sections  in  series  is  determined 
by  computer  and  compared  with  the  corresponding  poles  of  the  actual    line 
being  modeled.     For  either  an  open  or  short  circuited  termination,   the 
performance  specification  is  that  in  the  range  from  60  to  400  miles  long, 
the  percent  error  between  the  poles  of  the  model   and  those  of  the  actual 
line  shall   not  exceed  5  percent.     This  criterion  results  in  the  require- 
ment for  15  to  65  L-sections.     Thirty  sections  are  of  4  mile  "length" 
and  thirty  five  are  of  8  mile  "length". 

With  the  number  and  "length"  of  the  model   sections  determined, 
electrical   parameter  data  for  a  number  of  typical   extra  high  voltage 
(EHV)  lines  are  suitably  scaled  to  select  a  base  value  of  model   inductance 
and  capacitance  for  the  4  mile  and  8  mile  sections.     The  effects  of  non- 
transposition  are  modeled  by  incorporating  two  mutual    inductors  between 
the  phases.     Zero-sequence  inductance  is  modeled  by  a  specially  designed 
low  Q  inductor  whose  inductance  decreases  appropriately  with  increasing 
frequency. 

To  ensure  that  the  overall   quality  factor  of  the  model    remains 
greater  than  or  equal   to  that  of  actual   lines  for  all   frequencies  between 
60  Hz  and  3060  Hz;   low  loss,  powdered  molybdenum-nickel-iron  toroid  cores 
are  selected  for  trie  inductors,  and  low  dissipation  factor  polystyrene 
capacitors  form  the  basic  capacitance  of  each  section.     With  the  designed 
quality  factor  of  the  model   exceeding  that  of  most  EHV  lines  in  operation, 
the  Q  of  the  model    is  brought  into  correspondence  with  actual    lines    (at 
60  Hz)  by  suitably  inserting  resistance  into  each  section. 
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CHAPTER  1 
INTRODUCTION 

1.1  Transmission  System  Simulators 

The  three  phase  power  transmission  line  model  to  be  designed 
will  form  part  of  a  power  transmission  system  simulator  (TSS)  for  use 
in  the  M.I.T.  Electric  Power  Systems  Engineering  Laboratory.  Knowledge 
of  the  transient  behavior  of  power  system  networks  is  essential  to  the 
system  designer  and  planner  in  order  to  correctly  dimension  or  rate 
the  various  components  of  a  power  system.  Circuit  breakers  and 
insulators,  for  example,  cannot  be  correctly  sized  without  some 
knowledge  of  the  magnitude  and  duration  of  system  fault  currents  or 
switching  overvol tages.  With  765  kV  systems  already  in  operation, 
and  1100  kV  and  1500  kV  systems  proposed  for  the  future,  it  is 
particularly  important  to  know  the  magnitude  and  shape  of  over- 
voltage  transients  caused  by  routine  switching  operations.  Models 
for  1100-1500  kV  lines  still  in  the  design  stage  are  of  obvious  value, 
and  the  very  limited  access  to  extra  high  voltage  (EHV)  lines  in 
operation  makes  a  model  of  these  lines  highly  desirable. 

1.2  General  Model  Parameters  and  Constraints 

A  transmission  line  model  (t.l.m.)  is  only  one  part  of  a  TSS. 
Most  simulators  include,  in  addition,  models  of  power  system  generators, 
transformers,  compensating  reactors,  lightning  arresters,  circuit 
breakers,  ground  fault  neutralizes,  etc.  In  the  present  t.l.m. 
design,  the  requirement  for  interface  with  other  TSS  components  places 
certain  constraints  upon  the  model.  In  particular,  compatibility 
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with  existing  TSS  transformers  dictates  that  the  nominal  operating 
voltage  of  the  t.l.m.  be  138V  line-to-line  and  that  nominal  phase 
current  be  38.6  ma.  Other  required  features  for  the  t.l.m.  are  that 
it  be  capable  of  modeling  vertical,  horizontal,  transposed,  and  non- 
transposed  EHV  line  configurations  from  60  to  400  miles  in  length, 
have  quality  factors  exceeding  those  of  existing  EHV  lines  from  60 
to  3060  Hz,  model  the  variation  of  zero  sequence  inductance  with 
frequency  over  the  range  60  to  3060  Hz,  and  model  lines  in  the 
voltage  range  from  345  kV  to  1500  kV.  In  addition,  the  positive 
sequence  inductance  and  capacitance  of  the  t.l.m.  must  remain 
essentially  constant  over  the  range  60  Hz  to  3060  Hz,  and  inductance 
must  remain  constant  with  varying  current  in  the  range  from  0  to  at 
least  one  ampere-thus  ensuring  constant  model  inductance  up  to 
25  p.u.  fault  current. 
1.3  Features  and  Scope  of  Design 

One  of  the  more  unique  features  of  this  t.l.m.  design  is 
that  it  is  to  be  capable  of  modeling  non-transposed  power 

transmission  lines.  Numerous  transmission  line  models  in  simulators 

1-4 
are  presently  in  operation   ,  but  all  in  the  U.S.,  at  least,  are 

transposed  models.  A  literature  search  covering  the  period  subsequent 

to  1960  disclosed  only  a  limited  literature  on  t.l.m. 's  -  all 

relating  to  transposed  models.  Many  actual  lines  are,  in  fact,  not 

transposed.  Previously,  with  relatively  lower  voltage  transmission 

lines,  it  was  assumed  that  the  small  effects  of  non  transposition 

could  be  ignored.  Non- transposed  lines  were  mathematically  and 

physically  modeled  as  though  they  were  transposed,  and  model  results 
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were  considered  to  agree  satisfactorily  with  actual    lines.     Recent 

5 
studies     on  EHV  lines  indicate,  however,  that  the  effects  of  non- 
transposition  can  no  longer  be  ignored.     This  thesis  is  an  attempt 
to  design  a  t.l.m.  which  will   include  the  effects  of  non-transposition 
and  document  the  design  in  sufficient  detail   that  persons  other  than 
the  designer  can  construct  and  operate  the  model. 
1.4     Extensions  of  the  Design  Concept 

The  t.l.m.   design  discussed  in  this  thesis  is  intended 
specifically  for  the  modeling  of  electric  utility  power  transmission 

systems  for  which  the  X-/X,    ratio  (per  unit  length)  at  60  Hz  is  of 

2 
the  order  of  230,000-270,000  mi     (X~  =  positive-sequence  capacitive 

reactance;  X.    =  positive-sequence  inductive  reactance).     The  design 

concept,  however,   is  not  at  all   limited  to  utility  lines.     As  will 

become  clear  in  later  sections  of  the  discussion,  any  band  of 

Xc/X,    ratios  could  have  been  provided  -  had  it  been  desirable  or 

necessary.     With  the  proper  selection  of  model   components,  any  audio 

range  power  transmission  system  involving  inductive  and  capacitive 

coupling  could  be  accommodated.     Thus,  the  model   concept  is  equally 

suited  to  modeling,  for  example,  under  ground  power  cables,  long 

ship-board  or  industrial   plant  cable  runs  or  long  synchro  cable  runs, 

etc..     In  addition,  the  mathematical   approach  developed  in  Chapter  3 

for  the  solution  of  the  impedance  polynomial   has  application  in  any 

area  where  the  impedance  concept  has  validity  -  including,  for 

example,  the  mechanical   impedance  in  mechanical   vibrations  problems. 
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CHAPTER  2 
MODEL  CONFIGURATIONS 
2.1  General  Non-Transposed  Models 

A  section  of  a  general,   lumped -element,  three-phase 
transmission  line  is  shown  in  Figure  2.1.     The  electrical   para- 
meters of  tne  model   follow  directly  from  the  inductance, 
capacitance,  and  resistance  matrices  of  the  actual   line.     With  the 
matrices  of  the  line  defined  as 


[L]  = 


[C]  = 


aa   ab   ac 

Lab  Lbb  Lbc 

ac   be   cc 


C   -C  .  -C 
aa   ab   ac 

*Cab  Cbb  "Cbc 

•C   -C.   C 
ac   be   cc 


(2-1) 


and 


[R]  = 


Raa  Rab  Rac 
Rab  Rbb  Rbc 
Rac  Rbc  Rcc 


tne  correspondence  between  line  and  model  inductance  and  capacitance 
parameters  is  as  follows  (See  Appendix  A.l): 


Laa  =  Lll  +  L00 
Lbb  "  L22  +  L00 
Lcc  =:  L33  +  L00 


Cab  "  C12 
Cac  =  C13 
Cbc  ::  C23 


Lab  =  L12  +  L00 
Lbc  =  L23  +  L00 
Lac  "  L13  +  L00 


Caa  =  C12  +  C13  +  C10 

r   =  r       +  r   +  r 

ubb   L12   u23  u20 

Ccc  =  C13  +  C23  +  C30 


(2-2) 
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The  resistance  RQQ  in  the  model   represents  only  an  approximate  model - 
ing  of  line  inter-phase  resistances.     Rather  than  providing  three 
separate  inter-phase  resistance  circuits  in  the  model,  each  of  the 
three  line  inter-phase  resistances  is  approximated  by  the  average 
of  the  three  i.e., 


R        -   R        ^R        -    Rdb  +    RbC  +   RdC 
Kab  "'  Kbc  '  Kac  "  3 


■  V 


(2-3) 


So  that, 


CR]  = 


R 


R 


a  a   in 

Rm   R 
R_ 


m 
Rm 


(2-4) 


cc 

Comparing  line  equations  with  those  of  Appendix  (A.l),  it  is  found 
that: 


R00  =  Rni 


R„  =  Rn  +  Rm 

aa  11    m 

Rbb  =  R22  +  "Vi 

Rcc  =  R33  +  ^i 


(2-5) 


Tne  line  parameters  given  in  equation  (2-1)  include  the  effect 
of  ground  wires,  so  that  the  model  parameters  also  include  ground 
wire  effects.  The  mutual  inductances  shown  in  Figure  2.1  represent 
"ideal"  mutual  inductors  adding  no  self  inductance  to  the  lines. 
The  value  of  the  mutual  inductances  is  nominally  about  ten  percent 
or  less  of  the  self  inductances. 

The  circuit  of  Figure  2.1  can  be  transformed  by  simple 
algebra  (See  Appendix  A. 2)  to  the  configuration  of  Figure  2.2. 
The  circuit  of  Figure  2.2  is  as  general  as  that  of  Figure  2.1  but 
has  one  less  mutual  inductor.  In  Figure  2.2  the  transformed 
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quantities,  in  terms  of  Figure  2.1,  are: 

l'=l-L  L'  =L-L  f  =  r 

Lll   Lll   L13        L12   L12   L13        L10   L10 

L22   L22  L13         L23   L23   L13         L20   u20 

L33   L33  L13         L00    00   L13         L30   L30 

C12  =  C12  Rll  =  Rll  R00  =  R00 

C13  =  C13  R22  =  R22 

C23  =  C23  R33  +  R33  <A"4> 

If  the  line  being  modeled  is  of  horizontal  configuration;  then  for 

both  circuits: 

L12  =  (L12  "  L13^  =  L23  =  ^L23  "  L13* 

Hi  =  ^Lll  "  L13^  =  L33  =  ^L33  "  L13^ 
c  =  r   =  c  =  r 

L10   L10   u30   L30 

c    =  r  =  r'  =  r 

L12   L12   L23   L23 

Rll  =  Rll  =  R33  =  R33 
The  model  section  shown  in  Figure  2.2  (hereafter  referred  to  as 

"the  model")  is  the  configuration  selected  for  the  design.  While 

it  is  by  no  means  the  only  configuration  possible,  this  model  has 

a  number  of  advantages  over  other  possible  models  of  equal  generality. 

First,  it  should  be  pointed  out  that  in  order  to  economize  on  the 

number  of  cores  to  be  purchased,  the  mutual  inductors  are  wound  over 

tne  self  inductors  in  the  line.  The  configuration  of  Figure  2.2 

requires  four  cores,  and  incorporates  the  feature  that  no  self- 

inductor  will  be  required  to  physically  accept  more  than  one  mutual 

winding  (the  self  winding  leaves  sufficient  window  area  for  only  one 

mutual).  Second,  the  configuration  in  Figure  2.2  ensures  that  no 

phase  line  is  required  to  incorporate  the  resistance  of  more  than 
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one  mutual  winding.     (The  resistance  of  each  phase  line  must  be 

kept  to  a  minimum  to  keep  Q  high.)     Third,  when  modeling  horizontal 

lines,  the  mutual   inductances  of  the  Figure  2.2  configuration  are 

equal -thus  incorporating  into  the  model   the  symmetry  which  actually 

exists  in  the  line  being  modeled.     (This  would  not  be  true  of  some 

other  non-transposed  models.   ) 

2.2    Transposed  Models 

As  mentioned  in  Chapter  1,  it  is  common  practice  to  simplify 

the  analysis  of  non-transposed  lines  by  assuming  that  they  are 

transposed.     In  most  non-transposed  lines  the  values  of  the  three  self 

inductances  of  the  phase  conductors  are  approximately  the  same,  and  the 

values  of  all   inter-pnase  parameters  are  not  too  different  from  one 

another.     As  shown  in  Appendix  A. 3,  the  transformation  to  a  transposed 

view-point  "averages"  the  line  parameters  and  "re-distributes"  them 

uniformly  among  the  three  phases.     The  result  is  a  balanced,  symmetrical 

model  which  for  many  applications  is  a  good  approximation  to  what  is 

in  reality  a  non-transposed  line.     A  schematic  of  a  transposed  line 

model  section  resulting  from  the  transformation  of  Appendix  A. 3  is 

shown  in  Figure  2.3.     In  terms  of  the  general  model  of  Figure  2.1, 

the  parameters  of  Figure  2.3  are: 

CT=     C10  +  C20  +  C30  (A.5c) 


Cj=  °12  +  Cf  +  °13  (A-5d) 

R]-  R"  -  Rf  *  *33  (A-5e) 

l]=  L"  +  *f  *  '33  (A-5a) 

T  .  L12  +  L23  +  L!3  (A.5b) 
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2.2.1  Symmetrical  Component  Equivalent  Circuits  for  Transposed  Model 

In  Appendix  A. 3  it  is  shown  that  the  transposed  line  model 
can  be  reduced  to  two  symmetrical  component  models.  The  results  are 
repeated  here  for  continuity.  With  the  definitions 

RP.S.  =  RS        RZ.S.  =  RS  +  3R00       CP.S.  =  C0  +  3CM 

r     =  r^         I     =  I   -  I 
LZ.S.   ^0         LP.S.   LS   LM 

and  LZ.S.  =  LS  +  ZLl   +  3L00' 

the  transposed  line  equations  (A-7)  -  with  positive  sequence  currents 

only-reduce  to: 


9V 

a  -  R 

dz       Vs. 

h+l9.S. 

8Ia 

at 

b  =  R 
3z     P.S. 

!b  +  LP.S. 

aib 

at 

3V 

C  -  R 
3z    KP.S. 

"c+Lp.S. 

ai 
c 

at 

(A-8) 


And  The  Current  Law  equations  (A-ll)  reduce  to 
«...    3Va 


az    p.s.  at 

aib  avb 

az~  =  cp.s.  at"  ^A"12^ 

air  avr 
c  =  cr    c 


az    p.s.  at  * 

These  results  show  that  for  positive  (or  negative)  sequence  voltages 
and  currents,  the  transposed  line  may  be  viewed  as  three  identical 
single  phase  circuits  (the  120°  phase  difference  understood). 
Similar  conclusions  apply  to  zero  sequence  voltages  and  currents 
where  it  can  be  shown  (Appendix  A. 3)  that  under  these  conditions 
the  transposed  line  equations  reduce  to: 


and 
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3V  31 

3  =  R,  r     I     +  L- 


3z  Z.S.  a       UZ.S.   3t 

3V,  31. 

3^=RZ.S.  h  +  ll.S.JT  <A"9> 

3V  31 

'  =  R-7  c  I_  +  L- 


3z        z.s.  xc     uz.s.  at    • 


3Ia 

3z 

=  c     a 

UZ.S.  3t 

3Ib 

3z 

9Va 
=  C       a 

uz.s.  at 

3Ic 

3V 

=  C,  c      * 

(A-13) 


3z        z.s.  at 


Thus,  the  non-transposed  model   can  b6  approximated  by  a  transposed 
model,  and  the  transposed  model   can,  in  turn,  be  analyzed  into  two 
symmetrical   component  equivalent  circuits  -  one  for  positive-sequence 
and  another  for  zero-sequence. 
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CHAPTER  3 
NATURAL  FREQUENCY  ANALYSIS 
3.1  Methods  and  Criteria  for  Predicting  Model  and  Actual 
Line  Correspondence 

As  discussed  in  Chapter  2  and  Appendix  A,  the  transposed 
line  can  be  analyzed  into  two,  symmetrical  component,  single-phase 
equivalent  circuits.     A  single-phase  L-section  model   as  "seen"  by 
positive  (or  negative)  sequence  voltages  and  currents  is  shown  in 
Figure  3.1a.     A  similar  L-section  model   as  seen  by  zero-sequence 
quantities  is  shown  in  Figure  3.1b.     Figures  3.1a  and  3.1b  reflect 
the  basic  "L"  configuration  implied  by  the  three-phase  models  of 
Figures  2.1,  2.2  and  2.3,  but  "ir"  or  "T"  section  models  are  equally 
valid.     In  fact,  under  certain  circumstances  and  with  certain  loads, 
the  V  or  "T"  configurations  may  have  advantages  over  the  "L" 
models  -  these  circumstances  will   be  discussed  more  fully  in  later 
sections. 

One  of  the  purposes  for  which  the  transposed  line 
approximation  can  be  used  is  that  of  establishing  a  "math-model  - 
predictor"  for  comparison  of  an  actual   line  and  its  proposed  models. 
Some  of  the  early  questions  in  the  design  of  a  lumped  element 
t.l.m.   are:     How  many  sections  are  required?     How  "accurate"  is  the 
model  with  ten -twenty- forty  sections?     Should  the  sections  be  "L" 
or  "tt"?     What  should  be  the  "length"  of  each  model   section?     For 
what  frequency  band  must  the  model   "correspond"  to  the  actual   line? 
Exactly  what  is  it  that  should  "correspond"  between  model   and  line- 
to  what  "accuracy"?     A  method  of  theoretically  comparing  an  actual 
line  and  its  proposed  models  is  clearly  desirable.     If  a  t.l.m. 
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could  be  constructed  to  exactly  duplicate  in  magnitude  and  phase  the 
impedance  function  of  the  actual  line  for  all  frequencies,  then  the 
model  could  be  said  to  have  "perfect  accuracy".  Perfect  accuracy 
would  probably  require  more  lumped  element  sections  than  economically 
practical,  so  that  some  compromise  is  required. 

As  the  first,  and  perhaps  the  simplest,  approach  toward 
answering  the  questions  above,  this  design  uses  as  a  comparison 
"device"  the  poles  of  the  driving  point  impedance  functions.  The 
poles  of  a  given  length  of  actual  line  are  compared  to  the  poles  of 
a  proposed  m-section  model.  The  number  of  sections,  m,  in  the  model 
is  "adjusted"  until  the  frequency  difference  between  corresponding 
poles  of  the  model  and  actual  line  is  sufficiently  small.  In  each 
case,  the  actual  line  and  model  compared  are  the  transposed, 
symmetrical -component  equivalents  of  the  three-phase  systems.  This 
simplification  is  sufficient  for  design  analysis  since  it  is  not  the 
detailed  nature  (e.g.  "splitting")  of  the  poles  that  1s  sought,  but 
rather,  a  reasonably  accurate  estimate  of  their  mean  position  in  the 
frequency  spectrum.  The  mean  positions  can  then  be  compared  to 
estimate  the  difference  between  the  actual  line  and  model.  The 
transposed,  symmetrical -component  approximation  is  applied  to  both 
the  matnernatical  model  of  the  actual  line  and  that  of  the  t.l.m., 
so  that  errors  introduced  by  this  approximation  are  self  compensating 

In  analyzing  tne  driving  point  impedance  function  of  both 
the  line  and  models,  it  is  convenient  and  sufficient  to  consider 
two  possible  extremes  as  terminations.  Both  the  sources  and  loads 
connected  to  the  terminals  of  the  line  are  non-linear  devices  (at 
least  under  some  conditions  of  operation  or  fault)  so  that,  as 
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viewed  from  the  line,  the  terminals  of  the  line  may  simultaneously 
appear  as  a  high  impedance  to  some  transient  components  and  a  low 
impedance  to  other  components.  This  "high-low"  terminal  impedance 
characteristic  is  provided  for  in  the  analysis  to  follow  by  consider- 
ing both  open  and  short  circuit  terminations.  The  decision  as  to 
wnat  is  a  "sufficiently"  small  difference  between  the  line  and  the 
model  is  deferred  to  a  point  later  in  the  analysis  where  it  can  be 
made  to  appear  a  little  less  arbitrary. 
3.2  Driving  Point  Impedance  of  Continuous  Lines 

The  relations  describing  the  driving  point  impedance  of 
continuous  lines  are  well  known  and  will  not  be  derived  for  this 
discussion.  Readers  interested  in  the  details  of  the  derivation 
may  consult  any  good  text  on  transmission  lines  such  as  references 
7,  8,  or  9.  For  the  general  transmission  line  circuit  shown  in 
Figure  3.2,  the  pertinent  impedance  relations  for  a  length,  i9   are: 

Zs.c.  =  Zo  tanh  yl  (3_1) 


Mere 


Z    =  ln   coth  yl.  (3-2) 

O.C.     0 


Z    =  Impedance  at  Z  =  I   with  a  short 

5  •  w  • 


circuit  termi nation. 


Z    =  Impedance  at  Z  =  £  with  an  open 

1/  •  w  • 


circuit  termination 


Y  =     /(R  +  sL)sC  (3-4) 
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s  =  Complex  phase  constant 

R  =  Positive  or  zero-sequence  resistance  per  unit  length 

L  =  Positive  or  zero-sequence  inductance  per  unit  length 

C  =  Positive  or  zero-sequence  capacitance  per  unit  length 
Tne  shunt  conductance,  G,  has  been  ignored  in  these  relations.     In 
tne  short  circuit  case,  tne  poles  are  seen  to  occur  when 

cosh  y£  =  0.  (3-5) 

Setting  y  =  y    +  jy-,  equation  (3-5)  can  be  expanded  as  follows: 

cosh(yr  +  JY-jH  =  0 

cosh(Yr^)cosh(JYi5-)  +  sinh(Yr^)sinh(JY^)  =  0 

cosh(YJ<)cos  y-i  +  j  sinh(Y  £)sin  y.l  -  0. 
Thus,  for  cosh  yi  =  0  it  is  necessary  that 

COSh(Yr£)cos  y^SL  =  0,  (3-6) 

and 

sinh(Yr5<)sin  y^i  =  0.  (3-7) 

Case  No.   1 : 

Suppose  Yr  =  0 

Then  eq.    (3-7)   is  satisfied  and  eq.    (3-6) 

implies  yi  -  jy^l  =  +  j(2n+l)ir/2        n=0,l,2 (3-8) 

Case  No.  2: 

Suppose  Yr  t   0  but  y.t   =  +  ^-      n=0,l,2 

Then  eqs.  (3-6)  and  (3-7)  are  never  both  satisfied. 
Therefore  Case  No.  1  gives  the  only  solutions  for  y   which  result 
in  the  poles  of  equation  (3-1).  A  similar  development  apDlied  to 
equation  (3-2)  gives  for  the  open  line: 

yl  =   +  jniT      n  =  0,1,2 (3-9) 
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If  equations  (3-8)  and  (3-9)  are  expanded  and  solved  for  the  complex 
pnase  constant,  s,  there  results: 

S.C.  CL  t\-  ^C    LC 


so.c.  =  "   (fr)iV(fr)'-^  n-0.1,2....  (3-11) 

If  tne  line  shown  in  Figure  3.2  is  lossloss   (R  =  0)  then: 

Zs.c.  *  JZo  tan  &l 
Zo.c.  =  -  JZo  cot  » 

zo=yF 

3     =     go/QT    =  phase  constant 


ws, 

.c. 

.  +  (2n  +  1)ti 
—     2£  /LC 

r 

n-0,1,2.... 

Mo. 

.c. 

niT 
~  I  /LC 

n-0,1,2 

3, 

,3 

Dri 

vinq 

Point  Impedance 

Of 

L- 

-Section  Models 

(3-12) 

(3-13) 

(3-14) 

(3-15) 

(3-16) 

(3-17) 

A  schematic  diagram  of  an  L-section  model   for  modeling  the 
line  of  Figure  3.2  is  shown  in  Figure  3.3.     Both  figures  refer  to 
the  symmetrical   component  equivalent  circuits  of  their  respective 
three-phase  systems.     Thus,  both  figures  represent  a  "transposed- 
line-approximation"  to  their  respective  non-transposed,  three-phase 
counterparts.     In  the  case  that  the  L-section  line  is  shorted,  a 
clear  distinction  must  be  made  as  to  which  end  is  the  source  and 
which  is  the  termination.     In  the  first  case,  a  short  can  be 
placed  across  tne  last  capacitor-removing  it  from  the  circuit.     This 
case  is  referred  to  as   "inductor  input".     In  the  second  case,  the 
short  can  be  placed  to  connect  the  "last"  inductor  to  ground.     This 
case  is  referred  to  as  "capacitor  input". 
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3.3.1      Inductor  Input  L-Section  Mode]    (Shorted) 

The  inductor  input  case  is  that  shown  in  Figure  3.3.     The 
derivation  of  the  driving  point  impedance  expression  is  rather 
lengthy  and  is  shown  in  detail    in  Appendix  B.l.     The  lumped  element 
parameters  of  Figure  3.3  are  defined  in  terms  of  the  distributed 
parameters  of  Figure  3.2  as: 

Rs  =  St  (3-18) 

Ls  ■  £  (3-19) 

cs  ■  r-  •  <3-20' 

Where  "m"  is  the  number  of  sections  in  the  model.     The  result  from 

Appendix  B.l   for  an  m-section;   lossy  model   is: 

(Rs+Ls)sinh  mA 
Zs.c.   =  sinh  mA-sinh(m-l)A   '  (B-16) 

The  denominator  of  equation   (B-16)   is  a  polynomial   in  the  real 

variable  cosh  A  so  that  the  solutions,   [cosh  A]       . ,  to  the  equation 

root 

sinh  mA-sinh(m-l )A  =  0  yield  the  propagation  constants  in  the 

propagation  relation 

s2L  C  sRC 

[cosh  A]roQt  =  1   +  ^±    +    -f§.  .  (B-7) 

Solving  equation  (B-7)  for  s  gives: 


ss.c.  ■  -  (fr)  ±  V  (§/  -  fz^  [1"(cosh  X)™J         (B"26) 

If  the  model  of  Figure  3.3  is  lossless  (P.  =  0), 

juL  sinh  mA 
Zs.c.  =  sinh  mA-sinh(m-l)A    m=l »2,3 (B-17) 
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2i    r 
oj  L  C 

cosh  A  =  1 ^-  (B-18) 


ws.c.  =iyEc  ^-(cosh  AVoot]  •  m=1»2>3  <B-27> 

For  an  example  application,  suppose  their  are  seven  sections 

in  a  proposed  model.     Then  m=7;  there  are  six  values  of  [cosh  A] root» 

and  six  values  of  s  or  w         .     The  six  poles  of  the  model   given  by 

equations   (B-26)  or  (B-27)  are  compared  to  the  first  six  poles  of  the 

actual   line  given  by  equations   (3-10)  or  (3-16).     The  difference  in 

frequency  between  two  corresponding  poles  is  a  measure  of  the  error 

of  the  model   for  that  pole  number.     Equations   (B-16),   (B-26),  and 

(3-10)  were  implemented  on  the  computer,  as  shown  in  Appendices  C  and 

D,  for  various  numbers  of  sections  up  to  about  forty.     A  plot  of  the 

resulting  pole  error  vs.   pole  number  is  shown  in  Figure  3.4a.     The 

percent  error  in  Figure  3.4a  is  defined  as 

I   [equation   (B-26)-I  [equation   (3-10)] 
*E™--1°°- Inequation  (3-10)J •  <3-21> 

The  parameters  used  in  equations   (B-26)  and  (3-10)  for  Figure  3.4a 
are  the  positive  sequence  values  of  a  765  kV  "typical   aluminum 
tower"  configuration: 

I  =  100  miles  (3-22) 

R  =   .0251   ohms/mi  (3-23) 

C  =   .021  uf/mi  (3-24) 

L  =  1.45  mh/mi    .  (3-25) 

The  pole  number  shown  in  Figure  3.4a  is,  of  course,  a  discreet  integer 
variable  with  pole  no.   1   corresponding  to  the  pole  of  lowest 
frequency  in  either  the  actual   line  or  model.     The  curve  for  a  given 
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number  of  sections  in  the  model   is  a  series  of  discreet  points  but  is 
drawn  continuous  to  avoid  confusion.     Figure  3.4a  shows,  for  example, 
that  at  least  twelve  sections  are  required  to  model   the  line  if  it  is 
desired  to  keep  all   of  the  first  five  poles  of  the  model  within  five 
percent  error. 

If  the  typical   parameter  values  of  equations   (3-22)-(3-25) 
are  substituted  directly  into  equations   (3-10)  or  (B-26),  it  is  found 
that  the  term  jr  1S  of  order  10'  whereas  the  terms 

(2n+1)2Tr2     nvt  2m2       n   „.    ,,         -, 

T or  ~T~    L1-00^  ^root-1 

4JTLC  JTLC  r00t 

7  8 

are  of  order  10     or  10     for  even  the  lowest  values  of  m  or  n.     This 

suggests  that  for  purposes  of  a  pole  error  analysis,  the  resistance 

terms  can  be  neglected  when  considering  the  positive  sequence  of  the 

lines.     This,  in  turn  suggests  that  a  pole  error  analysis  using  the 

lossless  equations   (B- 17),   (B-27),  and  (3-16)  should  give  results 

very   similar  to  the  lossy  case  shown  in  Figure  3.4a.  The  results 

for  the  lossless  case  paralleling  the  analysis  leading  to  Figure 

3.4a  are  shown  in  Figure  3.4b.  The  percent  error  in  Figure  3.4b 

is  defined  as 


(3-26) 


As  expected,  Figures  3.4a  and  3.4b  are  practically  identical  so  that 
losses  will  be  ignored  throughout  the  remainder  of  the  pole  error 
analysis.  Note  that  equation  (3-26)  is  independent  of  £,,  L,  and  C, 
so  that  the  pole  error  curves  of  Figure  3.4b  apply  to  any  line  and 
model  provided  the  model  is  "inductor  input",  L-section,  and  both 


ijh  n-(cosh  x)root]  - 

(2n+l)7T 
21  fi£ 

%  Error  = 

100                           (2n+l)TT 
21  /LC 
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line  and  model  are  terminated  in  a  short.  The  pole  frequencies 
corresponding  to  the  first  fifteen  poles  of  a  shorted,  typical 
765  kV  line,  for  various  line  lengths  are  shown  in  Table  3.1. 

Equations  (B-26)  and  (3-10)  for  the  configuration  of  Figure 
3.3  were  also  run  using  zero  sequence  values  of  a  typical  100  mile 
765  kV  line.  The  pole  error  results  for  this  case  are  shown  in 
Figure  3.5  with 

I   =  100  miles  (3-27) 

R  =  .51  ohms/mile  (3-28) 

C  =  .0143  uf/mile  (3-29) 

L  =  4.2  mh/mile  .  (3-30) 

The  percent  error  in  Figure  3.5  is  as  defined  by  equation  (3-21).  A 
comparison  of  Figures  3.4a  and  3.5  shows  that  the  pole  error  for  zero 
sequence  parameters  is  practically  identical  to  the  error  for  positive 
sequence  parameters.  In  addition,  Figure  3.5  is  practically  identical 
with  the  lossless  results  shown  in  Figure  3.4b.  Thus,  throughout  the 
remainder  of  the  pole  error  analysis,  the  lossy  line  and  zero  sequence 
pole  error  will  be  assumed  to  be  the  same  as  the  lossless  line  error 
of  any  given  configuration. 
3.3.2  Capacitor  Input  L-Section  Model  (Shorted) 

If  the  model  of  Figure  3.3  is  turned  around  so  that  the 
capacitor  is  at  the  input,  the  result  (lossless)  is  shown  in  Figure 
3.6.  The  derivation  of  the  driving  point  impedance  of  this  line  is 
carried  out  in  Appendix  B.2.  The  analysis  considers  only  the  lossless 
case  since,  as  discussed  in  section  3.3.1,  the  effect  of  resistance 
can  be  neglected  in  the  pole  error  analysis.  The  results  from 
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Table  3.1 
Short  Circuit  Pole  Frequencies   (Pos.   Seq.) 


Pole 
No. 

60  Mi 

100  Mi 

200  Mi 

300  Mi 

400  Mi 

1 

760  Hz 

456  Hz 

228  Hz 

152  Hz 

114  Hz 

2 

2,280 

1,370 

685 

456 

343 

3 

3,800 

2,280 

1,140 

760 

570 

4 

5,345 

3,200 

1,600 

1,054 

800 

5 

6,840 

4,100 

2,050 

1,370 

1,025 

6 

8,335 

5,020 

2,560 

1,675 

1,257 

7 

9,850 

5,920 

2,960 

1,975 

1,480 

8 

11,390 

6,840 

3,420 

2,275 

1,710 

9 

12,900 

7,750 

3,880 

2,580 

1,940 

10 

14,400 

8,650 

4,325 

2,880 

2,160 

11 

15,950 

9,580 

4,790 

3,190 

2,390 

12 

17,410 

10,470 

5,235 

3,490 

2,620 

13 

18,950 

11,390 

5,695 

3,790 

2,840 

14 

20,500 

12,300 

6,150 

4,100 

3,075 

15 

22,000 

13,200 

6,600 

4,400 

3,300 

Table  of  Pole  Frequencies  for  Shorted,  Typical 
765  kV  Line  (Positive  Sequence-Lossless) 
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Appendix  B.2  are, 


jwL  sinh  mX 
s  

"s.c.  "  sinh(m+l  )X-sinh  mX 


Zc      ,         '     c*nh/m!lll_c4nh     ^  "1=1,2....  (B-39) 


w2L  C 
cosh  X  =  1 ^-^  (B-33) 

ws.c.   =l7rr^^cosh  AW   '       m=^2,3....  (B-40) 

Note  that  equation   (B-39)  is  not  the  same  as  equation   (B-17),  so  that 
the  values  of  (cosh  *)root  obtained  from  (B-39)  for  use  in   (B-40)  will 
be  different  from  the  roots  obtained  from  (B-17)  for  use  in  (B-27). 
Equations    (B-39),   (B-40),  and   (3-16)  were  solved  with  the  computer 
as  shown  in  Appendices  C  and  D.     A  plot  of  the  resulting  pole  error 
vs.   pole  number  is  shown  in  Figure  3.7.     Having  neglected  resistance, 
the  poles  of  the  model   are  given  by  equation   (B-40)  and  those  of 
the  actual   line  by  equation   (3-16).     The  percent  error  in  Figure  3.7 
is  defined  by  equation   (3-26).     Figure  3.7  shows  that  the  error 
behavior  for  the  capacitor  input  case  is  radically  different  from 
that  of  the  inductor  input.     The  capacitor  input  model   consistently 
under-estimates  the  actual   pole  frequency,  whereas  the  inductor  input 
model   over-estimates  the  frequency  of  the  lower  order  poles,  and, 
throughout  its  range,  yields  pole  frequencies  higher  than  the 
capacitor  input  model.     The  inductor  input  model   gives  higher  pole 
frequencies  due  to  the  fact  that  the  shorted  termination  removes 
the  last  capacitor  from  the  model -thus  reducing  overall  model   capacity 
and  increasing  the  model's  natural   frequencies.     Note  also,  that  the 
capacitor  input  model   errors  are  larger  at  lower  order  poles  and 
smaller  at  higher  order  poles  than  corresponding  lengths  of  inductor 
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input  model . 

3.3.3  L-Section  Model,  Open  Termination 

In  the  case  of  an  open  termination,  the  poles  are  unaffected 
by  which  terminals  are  considered  to  be  the  input  or  output,  and  only 
one  analysis  is  required.  The  schematic  of  the  open  circuit  case  is 
just  Figure  3.3  with  the  short  removed.  The  driving  point  impedance 
for  this  case  is  derived  in  Appendix  B.3.  The  results  are  (R  =  0): 


Zo.c.  =  juLs 


cosh  X  -   1  - 


sinh(m+l )X-sinh  mX 
2  sinh  mX-sinh(m+l )X-sinh(m-l )X 


"Vs 


oo 


o.c. 


ijh  [1"(cosh  A)root3 


(B-46) 
(B-18) 
(B-47) 


Where,  again,  use  has  been  made  of  the  fact  that, 

Ls  ■  r  •  cs  ■  r  •  Rs  ■  °. 

m  =  number  of  sections  in  the  model, 
I  =  length  of  the  line  being  modeled, 
and  L  and  C  are  the  electrical 
parameters  per  unit  length  of  any 
given  line. 
Equations   (B-46),    (B-47),  and  (3-17)  were  solved  with  the  aid  of  a 
computer  as  discussed  in  Appendices  C  and  D.     The  results  are 
shown  in  Figure  3.8  with  the  pole  error  defined  as 


%  Error  =  100 


m 

/^Cl-(coshX)root- 

nu 

mt 

(3-31) 


lAI 


As  in  equation  (3-26),  (3-31)  is  independent  of  I,   L,  and  C,  and 
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depends  only  upon  line  configuration.     Again,  the  values  of  the  roots 
are  different  from  what  they  were  for  either  of  the  preceeding  cases. 
The  pole  frequencies  corresponding  to  the  first  fifteen  poles  of  an 
open,  typical   765  kV  line,  for  various  line  lengths  is  shown  in 
Table  3.2. 
3.4     Driving  Point  Impedance  of  Tr-Section  Models 

Referring  to  the  line  shown  in  Figure  3.6,  it  is  seen  that 
the  L-section  model   can  be  converted  to  a  it  model  by  a  simple 
adjustment  of  the  first  and  last  sections.     If  the  capacitor  at 
the  mth  node  is  replaced  by  one  of  value  C  /2,  and  another  capacitor 
of  value  C  /2  is  added  at  the  "0"th  node,  the  result  is  a  Tr-section 
model.     A  similar  type  of  adjustment  could  be  made  using  inductors 
of  value  L  /2,  and  the  result  would  be  a  T-section  model.     The  it 
and  T-models  are  reciprocal   networks,  whereas  the  L-model   is  not. 
A  pole  error  analysis  for  TT-sections  follows  in  sections  3.4.1   and 
3.4.2.     The  analysis  for  T-sections  is  not  included  since  it  is 
assumed  that  it  will   add  little  to  all   the  cases  considered. 
3.4.1     TT-Section  Model    (Shorted) 

The  derivation  of  the  driving  point  impedance  relations  for 

the  shorted,  -rr-model   are  carried  out  in  Appendix  B.4  and  the  solution 

is  carried  out  in  Appendices  C  and  D.     The  result  for  the  shorted 

•rr-model   is  shown  in  Figure  3.9.     The  percent  error  in  Figure  3.9  is 

defined  as  in  equation  (3-26),  but  again  it  should  be  noted  that 

the  values  of  [cosh  A)       .    differ  for  each  line  configuration  and 

'root  3 

termination. 
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Table  3.2 


Open  Circuit  Pole  Frequencies 


Pole 
Uo. 

60  Mi 

100  Mi 

200  Hi 

300  Mi 

400  Mi 

1 

0  Hz 

0  Hz 

0  Hz 

0  Hz 

0  Hz 

2 

1,510 

905 

452 

315 

226 

3 

3,020 

1,810 

905 

604 

402 

4 

4,540 

2,720 

1,360 

906 

680 

5 

6,040 

3,620 

1,810 

1,205 

905 

6 

7,550 

4,540 

2,270 

1,510 

1,135 

7 

9,055 

5,440 

2,770 

1,810 

1,385 

8 

10,580 

6,350 

3,175 

2,120 

1,585 

9 

12,090 

7,250 

3,625 

2,420 

1,815 

10 

13,600 

8,175 

4,092 

2,720 

2,046 

11 

15,050 

9,050 

4,525 

3,070 

2,260 

12 

16,600 

9,998 

4,999 

3,325 

2,500 

13 

18,200 

10,900 

5,450 

3,640 

2,725 

14 

19,650 

11,800 

5,900 

3,940 

2,950 

15 

21,200 

12,700 

6,350 

4,240 

3,175 

Table  of  Pole  Frequencies  for  Open,  Typical 
765  kV  Line  (Positive  Sequence-Lossless) 
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3.4.2     u-Section  Model    (Open) 

The  open  circuit  case  is  derived  in  Appendix  B.5  and  solved 
in  Appendices   C  and  D.     The  results  are  as  shown  in  Figure  3.10.     For 
Figure  3.10,  the  percent  error  is  as  defined  in  equation   (3-31)  with 
due  recognition  that  the  values  of  (cosh  Mroot  in  this  case  are 
different  from  their  values  in  the  opened,  L-section  line. 
3.5     Conclusions 

1.  The  results  shown  in  Figures  3.4a,  3.4b,  and  3.5  indicate 
that  for  typical   EHV  lines  and  models,  the  percent  error  between  lines 
and  models  is  largely  independent  of  losses  in  either.     Thus,  the 
percent  error  for  either  positive  or  zero-sequence  is  the  same  for 

a  given  number  of  model  sections  and  a  given  type  of  termination. 

2.  The  results   indicated  by  Figures  3.4b  and  3.7  show  that 
due  to  the  non-reciprocal   character  of  the  L-section  model,  the  error 
between  model   and  line  depends  greatly  upon  which  terminals  of  the 
model   are  considered  input  or  output.     In  fact,  the  capacitor  input 
case  is  seen  to  be  considerably  inferior  to  the  inductor  input  case 
for  low  order  poles. 

3.  Comparison  of  the  shorted  L-model  with  the  shorted  tt- 
model   shows  that  the  L-model   has  a  rather  more  complicated  error 
pattern  than  the  ir-model .     The  L-model   first  over  estimates  the 
true  pole  position  for  low  order  poles  -  then  under  estimates  for 
higher  order.     The  ir-model   consistently  under  estimates  the  position 
of  the  actual   pole. 

4.  The  open  circuited  tt  and  L-models  give  practically 
identical   pole  error  plots.      (Actually,  each  Tr-model   has  one 
additional   pole  higher  than  the  highest  pole  of  the  corresponding 
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L-model . ) 

5.  Any  L-section  model   can  be  converted  to  either  a  tt  or  T 
model  by  a  simple  change  in  the  first  and  last  sections   to  obtain 
any  of  the  pole-error  patterns  of  preceeding  sections. 

6.  A  suggested  method  of  using  the  pole  error  plots  is 
shown  in  Figure  3.11.     Figure  3.11   shows,  as  an  example,  the  pole 
error  "curves"  of  a  twelve  section  L-model   as  taken  from  Figures  3.4b 
and  3.8.     A  five  percent  error  bound  is  suggested  by  the  dotted  lines 
at  the  top  and  bottom  of  the  figure.     The  region  between  the  solid 
curves  is  the  "error  band"  of  the  twelve  section  L-model.     If  the 
terminals  of  the  line  appear  as  a  predominantly  high  impedance  to 
frequencies  corresponding  to  the  first  five  poles,  the  line  error 

is  given   (approximately)  by  the  "open"  curve.     A  corresponding 
observation  applies  for  the  low  impedance  case.     The  dashed  curve 
mid-way  between  the  solid  curves  is  meant  to  indicate  the  "nominal 
error"  when  the  terminals  appear  as  neither  open  nor  shorted.     If 
twelve  sections  were  used  to  model   a  sixty  mile  line,  Figure  3.11 
shows  that  the  error  band  lies  within  the  five  percent  error 
bounds  up  to  about  the  fifth  pole.     As  seen  from  Tables  3.1   and  3.2 
the  fifth  pole  corresponds  to  about  6000  Hz  for  a  sixty  mile  line. 
One  could  say  then  -   "the  12-section  L-model   is  good  to  within  five 
percent  up  to  6000  Hz"   (for  a  765  kV  line).     An  error  plot 
similar  to  that  of  Figure  3.11   could  be  made  for  any  model   type 
composed  of  any  number  of  sections.     In  particular,  if  it  is 
desired  to  model   a  400  mile  line  to  "within  five  percent"  an 
inspection  of  Figures  3.4b,  and  3.8  shows  that  a  model  of  about 
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forty  sections  will  ensure  that  the  line  is   "good  to  within  five 
percent"  up  to  the  fifteenth  pole.     Tables  3.1   and  3.2  show  that 
the  fifteenth  pole  corresponds  to  a  frequency  of  about  3100  Hz 
(for  a  typical   765  kV  line).     Thus,  a  minimum  of  forty  ten-mile 
sections  would  be  required  to  model   a  400  mile  line  to  within  the 
five  percent  criteria. 

Tables  3.1   and  3.2  give  the  pole  frequencies,  specifically, 
for  a  typical   765  kV  line.     Actually,  as  will  be  shown  in  Chapter  4, 
corresponding  tables  for  most  other  EHV  lines  would  be  very  similar 
to  Tables  3.1   and  3.2   (within  about  ten  percent).     With  Tables  3.1 
and  3.2  as  representative  of  EHV  lines,  it  is  now  possible  to  make 
some  quantitative  decisions  concerning  the  number  of  sections  to  be 
built  for  the  t.l .m. 

7.  As  stated  earlier  in  this  report  -  as  a  minimum  -  the  t.l.m. 
is  to  model   EHV  lines  from  sixty  to  400  miles  in  length.     The  five 
percent  "specification"  referred  to  earlier  is  actually  a  result  of 
observations  such  as  those  in  Conclusion  6.     It  was  noted  earlier 
that  forty  sections   is  the  minimum  number  which  will   just  ensure 
that  a  400  mile  line-model   is  within  five  percent  to  3060  Hz.     This 
in  turn,  implies  at  least  forty  ten-mile  sections  for  the  model.     In 
addition,  it  is  proposed  to  model   all   line  lengths  between  sixty 
and  400  miles  by  changing  the  number  of  sections  in  the  model.     Thus, 
if  a  line  to  be  modeled  is  396  miles  long,  it  would  be  desirable  to 
set  up  a  model   "396  miles  long".     This  adjustment  could  be 
accomplished  with  the  ten-mile  sections  by  using  "taps"  to  adjust 
length.     However,  this  t.l.m.   design  must  model,  horizontal   or 
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vertical,  non- trans  posed  lines  -  this  too  requires  taps.     To  avoid 
having  to  tap  the  sections  for  length,  it  is  proposed  to  make  the 
sections   "fine"  enough   (i.e.  short  enough)  so  as  to  make  tapping 
unnecessary  for  length  adjustments.     Line  lengths  to  400  miles  can 
be  reached  to  within  one  percent  by  making  the  largest  sections 
have  only  eight-mile  length.     Thus,  the  396  mile  line,  e.g.,  could 
be  reached  (to  one  percent)  by  either  forty  nine  or  fifty  eight- 
mile  sections.     But  eight-mile  sections  are  not  "fine"  enough  to 
model   short  lines,  etc.,  etc.     Finally  then,  with  considerations 
such  as  those  above,  it  was  decided  that  the  t.l.m.  design  should 
consist  of  thirty  four-mile  sections  and  thirty-five  eight-mile 
sections.     This  selection  is  seen  to  be  quite  conservative  (by 

the  chosen  criteria)  and  provides  a  model  somewhat  "finer"  than 

1-4 
some  other  t.l .m. 's. 

8.  The  pole  error  analysis  is  seen  to  provide  a  means  of 
quantitatively  comparing  proposed  models,  and  is  useful  in  indicat- 
ing the  relative  merit  of,  say,  doubling  the  number  of  sections  - 
or  halving  them,  etc. 
3.6  Applications  and  Recommendations 

The  solution  technique  employed  to  solve  the  impedance 
functions  of  this  chapter  can  have  application  in  other  fields  of 
study.  The  denominator  of  each  of  the  impedance  functions  of  this 
chapter  corresponds  to  the  "frequency  polynomial"  of  mechanical 
vibrations  problems.  If  the  frequency  polynomial  is  set  equal  to 
zero,  it  becomes  the  "frequency  equation"  for  determining  the 
natural  frequencies  of  mechanical  systems.  Even  in  systems  with 
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only  two  degrees  of  freedom  the  frequency  equation  is  of  fourth  degree 
(quadratic  in  u  )  so  that,  in  continuous  systems  or  systems  having 
more  than  two  degrees  of  freedom,  complicated  solution  techniques 
must  be  employed  or  approximations  made.  The  technique  suggested 
in  this  chapter  and  employed  in  Appendices  C  and  D  is  a  means  of 
having  the  computer  generate  the  frequency  equation  and  then  solve 
it  for  the  natural  frequencies  of  the  system.  A  detailed  investiga- 
tion into  the  exact  class  of  problems  for  which  this  solution  technique 
may  have  validity  has  not  been  conducted  by  the  author,  but  it  is 
suggested  that  such  an  investigation  could  be  worthwhile. 
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CHAPTER  4 
PARAMETERS  OF  EHV  LINES  AND  MODELS 
4.1  Parameters  of  Representative  EHV  Lines 

In  order  to  determine  the  parameters  or  component  values  for 
a  t.l.m.,  it  is  first  necessary  to  know  the  parameters  of  typical 
lines  the  t.l.m.  would  be  expected  to  model.  The  t.l.m.  should  be 
capable  of  modeling  not  only  lines  already  in  existence  but,  in 
addition,  it  should  be  capable  of  modeling  proposed  and  even  hypothetical 
lines.  For  information  pertaining  to  the  parameters  of  existing  lines 
use  was  made  of  Report  No.  710,  Electrical  Characteristics  of  Trans- 
mission Lines,  of  the  American  Electric  Power  Service  Corporation 
(AEP)  .  Appendix  IV  of  Report  No.  710  describes  in  detail  thirty- 
four  different  EHV  transmission  line  configurations  in  use  by  the 
company.  From  the  tin' rty- four  configurations  a  sample  of  six  was 
selected  as  representative  of  types  of  configurations.  For  example, 
AEP  type  "H"  was  taken  as  one  representative  of  a  345  kV,  single 
conductor  per  phase,  single  ground  wire,  vertical  configuration. 
Type  "H-Special"  was  selected  since  it  is  a  common  line  and  represents 
a  form  of  triangular  configuration.  A  hypothetical  type  "39F"  was 
"constructed"  by  using  two  phase-conductors  per  bundle  spaced 
eighteen  inches  apart  -  this  was  taken  as  representative  of  a  two- 
conductor,  "vertical",  345  kV  line.  A  type  "41  A"  line  was  selected 
to  represent  a  345  kV,  two-conductor,  horizontal  configuration.  The 
AEP  "Typical  Aluminum  Tower"  configuration  was  taken  as  representative 
of  most  765  kV  lines  (horizontal).  A  hypothetical  765  kV  "H- 
Special"  was  constructed  by  scaling  the  345  kV  "H-Special"  conductor 
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distances  in  the  ratios  that  exist  between  the  type  "41  A"  and  the 
"Typical   Aluminum  Tower".     Finally,  two  hypothetical    1500  kV  con- 
figurations were  taken  from  reference  11   -  a  horizontal   configuration 
(Geometry  No.   2  and  Bundle  Drawing  No.   4),  and  a  triangular  line 
(Geometry  No.   3  and  Bundle  Drawing  No.   4). 

Typical   positive  and  zero-sequence  parameters  obtained  for 
these  representative  lines  are  shown  in  Tables  4.1   and  4.2.     The 
values  shown  in  Tables  4.1   and  4.2   (except  for  Xp/X.    and  Q)  were 

obtained  with  the  aid  of  a  computer  using  the  "Electric  Power  Systems 

1 2 
Engineering  Laboratory  (EPSEL)  Line  Parameters  Program"         This 

program  requires  as  input  the  geometrical  configuration  of  the 
conductors  of  a  given  line  (along  with  conductor  and  earth 
properties)  and  computes  the  inductance,  capacitance,  and  resistance 
matrices  of  the  actual  line  as  a  function  of  frequency.  In  addition, 
as  well  as  the  parameters  of  Tables  4.1  and  4.2,  the  program  computes 
and  prints  the  values  of  the  model  parameters  of  Figure  2.2.  The 
model  parameters  are  printed  on  the  basis  of  a  one-mile  long  section 
and  assumes  a  scale  factor  of  one  between  impedances  of  the  line  and 
model.  Model  parameter  values  at  sixty  Hz,  on  the  basis  described 
above,  for  the  eight  representative  line  models  are  shown  in  Table 
4.3  (earth  resistivity  is  taken  to  be  100  ohms-meters).  Referring 
to  Table  4.2  note  that,  except  for  the  1500  kV  lines,  the  positive- 
sequence  (or  average)  Q  at  60  Hz  is  about  twenty  or  lower.  The  high 
Q  values  for  the  1500  kV  lines  are  questionable  and  would  be 
difficult  to  duplicate  in  a  practical  model.  Consequently,  no 
attempt  will  be  made  to  even  approach  the  Q  values  suggested  by  the 
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345  kV 
H-Spec 


Type 
H 


765  kV 
H-Spec. 


Type 
41A 


1500  kV 
Horizontal 


1500  kV 
Triangular 


765  kV 
Horizontal 


Pos.  Seq.      Pos.  Seq.      Pos.  Seq. 


Table  of  Positive  Sequence  Parameters  for 
Some  Typical  EHV  Line  Configurations  (at  60  Hz) 

Table  4.1 


X, 


1/XC  XL  R  "L 

lor6  7-11   S^TTH  °-594^TTT  °-063  ^rfe  236 ,000   (mi .  )2 

39F  meg  ft  mile  mile  ' 

5.25  0.810  0.087  263,000 


5.65  0.751  0.087  236,000 


7.49  0.561  0.027  239,000 


6.82  0.627  0.054  234,000 


8.08  0.522  0.005  237,000 


8.41  0.494  0.005  240,000 


7.78  0.548  0.028  236,000 


58 


Zero  Seq.  Zero  Seq.  Zero.   Seq.  Pos.   Seq 

XL  R  1/XC  Q 


IT  1.59  ^r-  0.329  £r  3.92  -M-fr 

39F  mi  mi  meg  fi 


345  kV 
H-Spec, 


Type 
H 


765  kV 
H-Spec, 


Type 
41A 


1500  kV 
Horizontal 


1500  kV 
Triangular 


765  kV 
Horizontal 


1.85  0.395  3.63 


1.97  0.395  3.30 


1.66  0.418  4.65  21 


1.79  0.598  4.80  12 


1.10  0.131  5.83  104 


1.23  0.133  4.81  99 


1.64  0.522  5.31  20 


Table  of  Positive  Sequence  Q  and  Zero  Sequence 

Parameters  of  some  Typical   EHV  Line 

Configurations   (at  60  Hz). 

Table  4.2 
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Parameter 

Type 
39F 

LU 

1.83 

L22 

1.76 

1  ' 
L33 

1.59 

1  ■ 

L00 

0.736 

H2 

0.297 

1  ' 

L23 

0.147 

C 
L12 

3.6 

C 
°23 

3.5 

L13 

1.5 

C 

ho 

11.2 

c 

L20 

9.2 

c 

L30 

10.8 

Rll 

56.2 

R22 

60.4 

R32 

71.4 

K00 

89.0 

345  kV 
H-Spec. 

Type 

H 

2.21 

2.21 

2.16 

2.16 

2.05 

2.05 

0.928 

0.928 

0.022 

0.276 

-0.051 

0.168 

1.3 

2.5 

1.5 

2.6 

1.5 

1.2 

10.0 

9.5 

9.6 

8.0 

9.7 

8.9 

80.0 

80.0 

84.6 

84.6 

94.9 

94.9 

103. 

103. 

Table  of  Unsealed  Parameters  for  Models 

of  Representative  EHV  Lines.      (L1    in  mh/mi . , 

C    in  nf/mi . ,  and  R1    in  mft/mi . ) 

Table  4.3 


765  kV 
H-Spec. 

1.52 

1.49 

1.42 

0.975 

0.004 

-0.038 
2.2 
2.7 
2.7 

13.2 

12.2 

11.5 

16.5 

24.0 

40.4 
130. 


Parameter 

Type 

41A 

Lii 

1.81 

1  ■ 
L22 

1.78 

1  ' 
L33 

1.81 

1  ' 

L00 

0.886 

L12 

0.212 

1  ' 
L23 

0.212 

C 
L12 

2.4 

C 
°23 

2.4 

C 
L13 

0.7 

C 

13.2 

c 

L20 

11.9 

c 

u30 

13.2 

Ril 

51.5 

R22 

58.0 

R33 

51.5 

R00 

182. 
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1500  kV  1500  kV  765  kV 

Horizontal  Triangular  Horizontal 

1.52  1.30  1.60 

1.40  1.46  1.59 

1.52  1.30  1.60 

0.410  0.604  0.820 

0.147  0.066  0.214 

0.147  0.066  0.214 

2.7  3.1  2.95 

2.7  3.1  2.95 

0.7  3.3  0.76 

15.3  12.5  14.7 

14.8  13.3  13.1 

15.8  12.5  14.7 

6.2  4.7  26.0 

4.0  6.7  30.6 

6.2  4.7  26.0 

41.5  42.3  165. 


Table  of  Unsealed  Parameters  for  Models 
of  Representative  EHV  Lines.      (L'   in  mh/mi . , 
C    in  nf/mi . ,  and  R  in  mft/mi . ) 

Table  4.3  (Cont. ) 
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1500  kV  line  values.     With  a  Q  of  twenty,  the  765  kV  horizontal   line 
will  be  taken  as  the  minimum  standard  (in  Q)  to  be  achieved  by  this 
t.l.m.     This  means   (referring  to  Table  4.2)  that  each  of  the  sixty- 
five  sections  of  the  t.l.m.  must  have  a  positive  sequence  Q  of  at 
least  twenty  at  sixty  Hz.     The  positive  and  zero-sequence  Q  of 
transmission  lines  is  not  constant  with  frequency  and,  in  fact, 
increases  with  frequency.     Graphs  of  the  positive  and  zero-sequence 
Q  of  representative  lines  are  shown  in  Figures  4.1   and  4.2.     The 
positive-sequence  Q  is  seen  to  increase  steadily  in  the  range  to 
3060  Hz  while  the  zero-sequence  Q  becomes  nominally  constant  in  the 
range  from  1500-3060  Hz.     Each  model  section  must  provide  a  Q 
characteristic  equivalent  to  those  of  Figures  4.1   and  4.2  in  order 
to  accurately  model   the  actual   line.     If  the  model   section  is 
designed  to  meet  the  positive-sequence  Q  requirement  of  the  typical 
765  kV  line,  than  it  can  be  made  to  model   lines  of  lower  Q  by 
simply  inserting  resistance  in  the  phase  conductors. 

The  positive-sequence  inductance  of  EHV  transmission  lines 
is  practically  constant  over  the  frequency  range  of  interest,  but 
the  zero-sequence  inductance  is  not.     A  graph  of  the  variation  of 
zero  sequence  inductance  with  frequency  is  shown  in  Figure  4.3  for 
the  typical   765  kV  line.     The  zero-sequence  inductance  is  seen  to 
drop  markedly  with  increasing  frequency-parti cularly  in  the  low 
frequency  region.     Thus,  in  the  model   also,  the  zero-sequence 
inductance  must  correctly  decrease  with  increasing  frequency  and 
produce  (as  nearly  as  possible)  the  scaled  equivalent  of  Figure  4.3. 
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4.2.1  Line/Model  Scale  Factors  and  General  Model  Parameters 

Due  to  the  fact  that  the  t.l.m.  must  be  compatible  with  existing 
TSS  transformer  and  compensating  reactor  models,  the  base  voltage  of  the 
t.l.m.  must  be  138  V  line-to-line,  and  the  ratio  of  the  base  impedance 
of  the  model  to  that  of  a  typical  765  kV  line  must  be  0.355.  This  means 
that  the  base  current  and  power  of  the  t.l.m.  are  38.6  ma  and  9.2  watts 
respectively.  The  base  impedance  of  a  typical  765  kV  line  on  a  100  MVA 
base  is  5850  ohms.  The  base  impedance  of  the  model  is  2078  ohms.  With 
the  above  impedance  ratio  constraint  between  the  model  and  765  kV  line, 

B  model  _  355 


and 


LB 

actual 

RB 

model 

RB 

actual 

CB 

model 

=  .355, 


=  2.82  . 


CB  actual   *355 
To  model  the  typical  765  kV  line,  all  the  765  kV  parameters  of  Table 

4.3  must  be  appropriately  scaled  by  the  ratios  above.  For  example, 

Lio  of  Table  4.3  becomes  .564  mh/mi ,  or  2.26  mh/4  mi,  or  4.52  mh/8  mi. 

Thus,  the  4-mile  section  would  be  constructed  with  a  basic  2.26  mh 

inductor.     An  identical   procedure  applied  to  the  345  kV  type  "41  A" 

line  gives   (on  a  100  MVA  base): 

(ZB)345     =  U90O 
ZBm         _  2078  =  ,   ,. 


iy^    "  1190 

LAp  =3.1  mh/mi,  or  12.4  mh/4  mi,  etc. 
To  model   the  765  kV  line  the  inductor  in  the  four  mile  section  must 
be  2.26  mh  and  to  model   the  345  kV  line  it  must  be  12.4  mh.     This 
implies  a  very  wide  range  of  taps  on  the  model   inductors.     A  wide 
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range  of  taps  on  an  inductor  is  not  only  difficult  and  expensive  to 
construct,  but  has  the  further  disadvantage  that  the  Q  of  the  inductor 
gets  "tapped"  as  the  inductance.  Thus,  tapping  the  inductance  down  to 
one-sixth  of  the  total  inductance  cuts  the  Q  down  to  about  one-sixth 
of  what  it  was  for  the  total  inductor. 

A  procedure  for  solving  this  problem  is  to  choose,  for  lines 
other  than  the  typical  765  kV,  a  power  base  other  than  the  commonly 
used  100  MVA.  Again  suppose  that  the  345  kV  type  "41  A"  line  is  to  be 
modeled.  The  first  step  is  to  establish  an  impedance  ratio  between 
the  actual  line  and  model  such  that  L^p  for  the  345  kV  model  is  the 
same  as  that  already  provided  for  the  typical  (horizontal)  765  kV 
line.  In  other  words,  using  Table  4.3,  seek  a  scale  factor  F, 

F  -  gBJ  ™>del  ,  (4-D 

UB;  line 
such  that 


So  that 


and 


F^L22^345  =  ^L22^  765  hor.  *  ^4"2^ 

c       ^L22^765  hor.  .  .564    ,1C  ,A   ^ 

F"^^ OS"'316'  <«> 

£S.  .316  . 

V2 
ZB*  "  T1   •  (4'4) 

(.316)  V2 

Pro  =  7 —  =  18.1  MVA  .  (4-5) 

U  ZBm 

Thus,  when  modeling  a  345  kV  type  "41A"  line,  it  is  necessary  to  use 


But 


so  that 
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a  power  base  in  the  line  of  18.1  MVA  instead  of  the  usual   100  MVA. 
The  resulting  parameters  required  of  a  four-mile  model   section  are 
shown  in  Tables  4.4a  and  4.4b  respectively.     Note  in  Table  4.4a  that 
the  inductors  corresponding  to  phases   "a"  and  "c"  are  equal   -  reflect- 
ing the  symmetry  of  the  horizontal   lines.     This  is  not  the  case  for 
the  "vertical"  lines  of  Table  4.4b.     With  the  inductances  of  the  765 

kV  horizontal  model   as  a  base,  note  that  all   the  L   's  of  Table  4.4b 

a 

are  larger  than  2.28  and  all  the  L  's  are  smaller.  In  fact,  all  the 

L  's  of  Table  4.4b  are  larger  than  2.28  by  the  same  amount  (A)  that 
a 

all  the  L  's  of  Table  4.4b  are  smaller  than  2.28.  The  scale  factor 
to  achieve  this  "balance"  for  vertical  configurations  is  not  given 
by  equation  (4-3)  but  is  determined  as  follows: 
Using  Table  4.3,  seek  a  scale  factor  F  such  that 

^HAert.  "  ^L11^765  hor.  =  A 

F(L33)vert.  "  (L33*765  hor.  =  "  A* 
Solving  for  F  [using  the  fact  that  (L^)^  hor>=  (L^)765  hor<], 

f  -  2(L11)765  hor.  ,,  ,x 

r  ~  7P — rr-\ •  v4-bj 

^Lll'L33;vert. 
The  use  of  equation  (4-6)  ensures  that  the  unbalance  inherent  in 
vertical  lines  is  shared  equally  by  both  phases  "a"  and  "b"  in  the 
model.  A  "balanced  sharing"  is  desirable  since  a  factor  determined 
as  in  equation  (4-3)  gives  the  result  that  the  phase  "c"  model 
inductor  would  be  required  to  span  a  larger  range  than  the  phase  "a" 
inductor  -  thus  requiring  more  taps.  Use  of  equation  (4-6)  ensures 
that  the  phase  "a"  and  "c"  inductors  will  have  the  same  number  of 
taps.  Referring  again  to  Table  4.4a,  note  that  L,  for  the  two 


Parameter 

Type 

41A 

La 

2.29 

Lb 

2.26 

Lc 

2.29 

Ln 

1.12 

Lab 

0.27 

Lbc 

0.27 

Cab 

30.4 

Cbc 

30.4 

Cac 

8.8 

Can 

166 

Cbn 

149 

Ccn 

166 

Ra 

65 

Rb 

73 

Rc 

65 

R 
n 

230 

68 


1500  kV  1500  kV  765  kV 

Horizontal  Triangular  Horizontal 


2.34  2.16  2.28 

2.16  2.41  2.26 

2.34  2.16  2.28 

0.630  1.0  1.16 

0.23  0.11  0.30 

0.23  0.11  0.30 

28.1  30  33.2 

28.1  30  33.2 

7.0  32  8.6 

164  121  166 

154  129  148 

164  121  166 

10  7.85  37 

6.5  10.8  43 

10  7.85  37 

67  70  235 


Table  of  Scaled  Parameters   Implied  for  a 
4-Mile  Hodel   Section.     (L     in  mh/4  mi., 
C     in  nf/4  mi.,   R  in  mS2/4  mi.) 

Table  4.4a 
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Parameter 

Type 
39 

345  kV 
H-Spec. 

Type 
H 

765  kV 
H-Spec. 

La 

2.44 

2.36 

2.36 

2.34 

Lb 

2.34 

2.31 

2.31 

2.29 

Lc 

2.12 

2.20 

2.20 

2.2 

Ln 

0.98 

0.99 

0.99 

1.50 

Lab 

0.40 

0.24 

0.30 

0.006 

Lbc 

0.20 

-0.05 

0.18 

-0.06 

Cab 

43.2 

19.4 

37.3 

22.8 

Cbc 

42.0 

22.4 

38.8 

28.1 

Cac 

18.0 

22.4 

17.9 

28.1 

Can 

135 

149 

142 

137 

Cbm 

111 

143 

119 

127 

Ccn 

129 

145 

131 

119 

Ra 

75 

86 

86 

20 

Rb 

80 

91 

91 

37 

Rc 

95 

102 

102 

62 

R 
n 

117 

110 

no 

201 

Table  of  Scaled  Parameters   Implied  for  a  4-Mile 
Model   Section.      (L     in  mh/4  mi . ,  C     in  nf/4  mi . ,   R     in  mft/4  mi.) 

Table  4.4b 
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1500  kV  lines  is  not  2.26  mh  as  would  be  found  had  equation  (4-3)  been 

used  to  determine  the  scale  factor.  In  fact,  still  another  technique 

has  been  employed  to  determine  the  scale  factor  of  these  two  lines. 

This  third  technique  will  be  discussed  later  in  this  chapter. 

Some  remarks  are  in  order  at  this  point  concerning  the  model 

zero-sequence  inductance,  l_z  s  ,  and  the  inductor  L  of  Table  4.4a. 

To  obtain  the  zero-sequence  characteristic  required  of  a  4-mile  model 

section  for  modeling  a  765  kV  line,  for  example,  the  graph  of  Figure 

4.3  is  simply  scaled  by  the  factor  (0.355)(4.0)  =  1.42;  the  result 

is  shown  in  Figure  4.4.  The  model  zero-sequence  inductance  is  required 

to  fall  from  6.2  mh  at  60  Hz  to  about  4.0  mh  at  1000  Hz.  The  point 

to  note  is  that  L  of  Tables  4.4a-b  is  not  the  zero-sequence  inductance 

of  the  model  per  se,  but  L  is  the  component  which  must  be  designed 

to  produce  the  characteristic  shown  in  Figure  4.4  if  a  transposed 

765  kV  horizontal  line  is  to  be  modeled.  The  detailed  design  of  a 

suitable  L  for  this  t.l.m.  is  being  carried  out  by  Mr.  William  Feero 
n 

as  part  of  his  thesis  to  be  completed  in  September,  1971.  Thus,  the 
design  of  L  is  being  carried  out  in  parallel  with  (but  slightly 
behind)  the  design  of  the  remainder  of  the  t.l.m. 

In  preparation  for  the  design  of  the  model  inductors, 
excluding  consideration  of  the  1500  kV  lines  of  Tables  4.4a  and  4.4b, 
it  was  decided  that  the  inductance  values  shown  in  Table  4.5  should 
be  provided  in  the  4-mile  sections  in  order  to  model  most  EHV 
lines. 
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mh 
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,08 

mh 

2. 

,21 

2. 

,13 

2. 
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2. 
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2, 

,31 

2. 

,23 

2, 
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2. 

,28 

2. 

.41 

2. 

,33 
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Table  4.5 
Self  Inductance  Designed  Into  the  4-Mile  Sections 

2.24  mh 

2.29 

2.34 

2.39 

2.44 

2.49 

Early  in  the  design  of  the  t.l.m.   the  1500  kV  lines  were 
generally  ignored  as  beyond  the  modeling  capabilities  attainable  by 
this  design.     With  respect  to  the  quality  factor,  the  preceeding 
statement  is,  indeed,  true.     However,  in  all  other  respects  it  is 
possible  to  model   these  lines  also-with  a  proper  choice  of  scale 
factor.     Note  in  Table  4.5,  under  L,  the  value  2.26  mh  originally 
provided  for  modeling  the  typical   765  kV  line.     Equation  (4-3)  shows 
the  Lpo  (unsealed  L  )  of  a  345  kV  line  being  scaled  to  the  2.26  mh 
inductor  previously  provided.     The  2.26  mh  value  is  not  unique,  and 
any  inductance  parameter  of  a  line  to  be  modeled  could  be  scaled 
to  any  one  of  the  inductances  provided  in  Table  4.5.     In  fact,  to 
properly  fit  the  two  1500  kV  lines,  the  L    of  the  horizontal   1500 
kV  line  was  scaled  to  2.16  mh  and  the  triangular  to  2.41  mh. 
(The  extremes  of  the  existing  L    were  used  in  order  to  keep  the 

resulting  values  of  (Oicqq  and  (La^l500  witn^n  "reacn"  of  tne 
ranges  provided  in  Table  4.5.)  The  above  procedure  resulted  in 
an  impedance  scale  ratio  of  0.414  for  the  triangular  line  and  0.384 
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for  the  horizontal.  The  resulting  power  bases  to  be  applied  to  the 
actual  lines  are  448  MVA  for  the  triangular,  and  417  MVA  for  the 
horizontal.  Thus,  the  three  scaling  examples  provided  in  this 
section  illustrate  the  considerable  flexibility  of  the  "Floating 
Power  Base"  technique.  The  values  of  self  inductance  to  be  designed 
into  the  8-mile  sections  are  obtained  by  simply  doubling  the  values 
in  Table  4.5. 
4.2.2  Selection  of  Model  Mutual  Inductance  Parameters 

Tables  4.4a  and  4.4b  indicate  that  typical  values  of  L  . 
range  from  0.006  mh/4  mi  to  0.4  mh/4  mi,  and  those  for  L.   from 
-0.06  mh/4  mi  to  0.27  mh/4  mi.  Providing  this  wide  range  of  values 
is  neither  practical  nor  necessary  for  the  following  reasons: 

1.  To  cover  the  range  implied  above  would  require  too 
many  taps  on  the  mutual  windings  and  too  many  values 
brought  to  the  front  panel  of  the  section.  As  few 

as  six  "tap  values"  per  component  in  Table  4.4  means 
that  120  lines  must  be  brought  to  the  front  panel. 
(It  is  proposed  to  limit  the  number  of  taps  on  any 
component  to  six.) 

2.  In  those  cases  where  the  mutual  inductance  value 
is  1/20  or  less  of  the  self  inductance  upon  which  it 
is  wound,  it  would  be  nearly  impossible  to  measure 
the  effect  of  the  mutual  inductance.  Therefore, 
those  values  of  mutual  inductance  amounting  to  1/20 
or  less  of  the  self  inductance  in  the  same  line  will 
simply  be  ignored.  Note  that  this  rule  leaves  the 
negative  mutuals  ignored. 
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3.  Discounting  the  extremes  mentioned  above,  most 
mutual   values  appear  to  be  in  the  range  of  from 
0.17  mh  to  about  0.27  mh  -  the  single  exception 
being  the  far  removed  0.40  mh  of  the  hypothetical 
type  "39F"  line.     Mutual   taps  actually  provided  in 
the  model  section  should  reflect  this   "crowding"  in 
the  mid-band  by  the  provision  of  relatively  more 
taps  in  this  band. 

4.  Horizontal   lines  are  common,  and  it  is  seen  that 
for  horizontal   lines  the  mutual   inductances  are  equal. 
The  selection  of  taps  to  be  provided  in  the  model 
section  should  preserve  this  balance  if  possible. 

Based  upon  the  arguments  above,  the  values  of  mutual   inductance 
to  be  designed  into  the  4-mile  model   sections  are  shown  in  Table 
4.6  -  those  for  the  8-mile  sections  are  shown  in  Table  5.5. 


Table  4.6 

Lab(mh) 

Lbc(mh) 

.15 

.15 

.18 

.18 

.20 

.20 

.23 

.23 

.25 

.25 

.33 

.32 

Table  of  Mutual   Inductance  Values 
Designed  into  the  4-Mile  Sections 
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4.2.3  Selection  of  Model  Resistance  Parameters 

The  resistance  ranges  implied  for  a  4-mile  model  section  are 
shown  in  Tables  4.4a  and  4.4b  in  which  the  figures  represent  the 
total  resistance  required  of  a  section  from  input  terminal  to  output 
terminal.  Thus,  the  resistance  indicated  includes  the  resistance  of 
the  self  and  mutual  inductance  windings  discussed  above.  The  amount 
provided  for  operator  adjustment  is  the  difference  between  winding 
and  connection  resistance  and  the  values  indicated  in  the  tables. 
The  resistance  of  the  self  and  mutual  windings  must  obviously  be  less 
than  the  resistance  values  in  Tables  4.4  if  the  4-mile  model  section 
is  to  provide  a  correct  modeling  of  line  Q.  The  exact  value  of  the 
resistances  to  be  designed  into  the  4-mile  section  cannot  be  deter- 
mined until  the  resistance  of  the  self  and  mutual  windings  is  known. 
The  setting  of  the  operator  adjusted  resistance  parameters  to  be 
provided  is  deferred  to  Chapter  7. 

4.2.4  Selection  of  Model  Capacitance  Parameters 

The  capacitance  parameters  implied  for  a  4-mile  section  are 
also  indicated  in  Tables  4.4a  and  4.4b.  With  a  wide  range  of 
capacitor  sizes  readily  available  from  manufacturers,  a  bank  of 
six  capacitors  will  be  provided  for  each  capacitance  parameter  for 
a  total  of  36  capacitors  in  each  section.  The  capacitors  of  each 
bank  are  to  be  sized  so  as  to  span  the  ranges  indicated  with  "room" 
on  either  side.  With  six  capacitors  per  capacitance  parameter, 
about  4!  values  can  be  presented  to  cover  the  range  of  the  parameter. 
Tables  of  the  actual  capacitor  banks  purchased  for  each  parameter 
are  presented  in  Chapter  6. 
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4.3     Summary  Remarks 

1.  The  parameters  of  the  six  non-1500  kV  lines  indicated 
in  Tables  4.1,  4.2,  and  4.3  are  taken  to  "span"  the  range 
of  typical   EHV  lines   (including  500  kV  lines). 

2.  No  attempt  is  made  to  model  positive  sequence  line 
Q's  greater  than  20. 

3.  Six  model  values  spanning  the  range  of  EHV  line 
configuration  "types"  are  provided  for  each  parameter. 

4.  Only  three  scaling  techniques  were  illustrated,  but 
any  line  to  be  modeled  could  be  scaled  to  any  one  of  the 
existing  model  values  in  an  effort  to  make  diverse  lines 
fit  the  model . 

5.  A  scale  factor,  F,  obtained  for  inductance  and 
resistance  scaling  is  used  as  1/F  for  capacitance 
scaling. 

6.  Values  of  mutual  inductance  less  than  five  percent 
of  the  self  inductance  in  the  same  line  are  ignored. 

7.  Parameters  indicated  in  Tables  4.4a  and  4.4b  for  a 
4-mile  model  section  are  simply  doubled  to  obtain  the 
values  for  an  8-mile  section. 

8.  The  details  of  the  return  path  inductor  design,  L  , 
will  appear  in  a  subsequent  report  (by  a  different  author) 
due  in  September,  1971. 
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CHAPTER  5 
INDUCTOR  DESIGN 
b.l  Model  Inductor  Requirements 

The  model  inductance  values  required  of  a  4-mile  section  are 
shown  in  Tables  4.5  and  4.6.  The  inductance  values  for  an  eight-mile 
section  are  simply  double  those  of  Tables  4.5  and  4.6.  In  addition 
to  simply  providing  the  inductance,  the  inductors  in  the  model  sections 
must  have  Q  considerably  greater  than  the  Q  characteristic  for  the 
"best"  line  to  be  modeled.  The  Q  characteristic  of  the  765  kV 
Horizontal  line  shown  in  Figure  4.1  is  taken  as  the  minimum  standard 
to  be  met  by  the  model  sections.  The  Q  characteristics  of  the 
individual  model  inductances  must  exceed  the  requirements  of  the  over- 
all section.  While  the  line  and  section  require  a  Q  of  twenty  at 
60  Hz,  the  inductors  should  have  a  Q  of  perhaps  40  at  60Hz  -  and  a 
similar  excess  over  the  entire  frequency  range.  The  characteristic 
shown  in  Figure  4.1  is  the  "average"  Q  of  the  three  phases  so  that 
it  is  not  the  exact  characteristic  to  be  tracked  by  any  of  the 
three  self  inductors  of  the  model  even  if  this  tracking  were  possible. 
Tne  solution  of  the  "Q  Tracking  Problem"  in  this  design  is  to  design  or 
select  an  inductor  whose  Q  characteristic  has  the  general  shape  shown 
in  Figure  4.1  but  whose  Q  value  is  approximately  fifteen  to  twentv- 
five  points  higher  along  the  range.  Any  excess  Q  of  the  model  can 
be  handled  by  simply  inserting  an  additional  small  resistance  in 
series  with  the  model  inductors. 

In  addition  to  the  Q  characteristic  per  se,  the  model 
inductors  must  produce  the  required  Q  at  the  yery   low  values  of  base 
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current  selected  for  the  model,  and  must  have  an  inductance  that  remains 
essentially  constant  up  to  at  least  25  p.u.   fault  current.     The 
inductance  of  all   the  phase  and  mutual   inductors  must  also  remain 
constant  with  frequency  -  only  the  return  path  inductor,  L   ,  is  to 
vary  with  frequency  to  provide  the  zero  sequence  inductance  characteristic, 
5.2     Inductor  Core  Selection  and  Core  Properties 

Tne  inductor  characteristics  described  in  section  5.1  are 
determined,  to  a  very  large  extent,  by  the  core  material  of  the  inductor. 
Tne  tneory  leading  to  tne  final   relations  describing  inductor  behavior 
is  rather  involved  and  will   not  be  included  in  this  report  -  only  the 
more  significant  relations  will   oe  presented  to  indicate  the  overall 
procedure.     Readers  interested  in  a  more  complete  background  leading  to 
tne  relations  presented  herein  should  consult  references  13-19. 

A  number  of  core  materials  and  configurations  were  investigated 
to  meet  the  specifications  described  above,  but  most  were  readily  found 
to  be  unsuitable  and  rejected  early  in  the  design.     A  potted  ferrite 
core,   for  example,  meeting  all  other  requirements,  was  rejected  as  too 
expensive  (Ferroxcube  series  6656).     Two  alternatives  were  considered 
in  some  detail   -  oriented  silicon  steel   "C"  cores  and  molybdenum- 
nickel-iron  powdered  toroids.     The  oriented  silicon  was  rejected  when 
tne  model  base  current  was  finally  set  at  the  relatively  low  value  of 
38.6  ma.     At  this  low  current  the  oriented  silicon  could  not  produce 
tne  desired  Q  characteristic  at  a  sufficiently  low  cost  or  weight. 
Basically,  the  Q  of  oriented  Fe-Si   inductors  increases  with  increasing 
current  (up  to  saturation),  whereas  with  molybdenum-nickel,  Q 
decreases  with  increasing  current.     Thus,  the  Mo-Ni   cores  nave  their 
highest  Q's  at  low  currents  and  the  Fe-Si   cores  have  their  lowest  Q's 
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at  low  currents.     The  basic  relation  governing  the  Q  of  an  inductor  can 
be  written  as   (for  the  frequencies  of  interest  here) 

Q  =  — 41 o-  (5-]> 

r  .  r  +  ri* 

dc  ac 

Where 

L  =  Inductance  in  henries 

co  =  2ir  x  frequency 

R.     =  dc  resistance  of  the  winding   (ft) 

Ra„  =  equivalent  core  loss  resistance  (ft) 
ac 

2 

With  these  definitions  it  is  seen  that  R  I  corresponds  to  both  the 

ac  ~— ~ ~~ 

hysteresis  and  eddy  current  losses  (See  reference  17  p.p.  58ff . )  and 

thus  represents  the  total  core  loss.  Ra  .  the  core  loss  coefficient, 

ac 

18 
is  itself  a  function  of  current.     For  powdered  cores  the  manufacturer 

provides  data  on  R      explicitly,  whereas  for  Fe-Si   the  entire  term 

ac 

R     I     is  provided.     If  the  data  provided  by  the  manufacturer  is 
ac 

converted  to  the  terms  of  equation  (5-1),  it  is  found  that  for  Fe-Si 
the  core  loss,  W,   is  given  by  (approximately) 

W  =  k^uOl1-7,  (5-2) 

whereas  for  the  powdered  Mo-Ni   the  loss  is  given  by   (approximately) 

W  =  RacI2  =  k2(w)I2'4  .  (5-3) 

In  equations   (5-2)  and  (5-3),  k,(to)  and  k2(w)  are  core  loss  constants 
and  are  a  function  of  frequency.     If  equations   (5-2)  and  (5-3)  are 
suustituted  in  equation  (5-1)  there  results,  respectively, 

0  = S»y (5-4) 

gFe-Si       „     T2   .    .    f    u1.7     •  ^D  *> 

RdcI     +  k^oojl 

and 
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QMo-Ni  =  I     T2  "  .    J    H2.4  •  (5_5) 

dc  ty   ' 

Equations   (5-4)  and  (5-5)  apply  only  in  the  region  below  saturation 

and  suggest  the  reason  why  the  Q  of  the  Fe-Si  core  increases  with 

current,  whereas  in  the  Mo-Ni   core  Q  decreases.     The  Q  calculations 

were  conducted  for  a  number  of  oriented  Fe-Si   cores   (Arnold-AH-393, 

AH-34,  AL-59,  AH-95,  AH-86,  AH-31),  maximized  with  respect  to  air  gap, 

with  the  result  that  cores  up  to  eight  pounds  each,  costing  up  to 

$39.00  each,  could  not  produce  the  required  Q  characteristic  at  the 

base  current  levels  required. 

Several  Mo-Ni   cores  could  meet  the  characteristic,  but  the 

smallest  (in  terms  of  those  offered  by  the  Arnold  Engineering  Co.  or 

i-lagnetics  Corp.)  and  least  expensive  unit  which  could  be  made  to 

produce  the  desired  characteristic  with  a  Q  of  40  at  60  Hz  was  Arnold 

part  no.  A-218259-2.     Tiie  Q  characteristic  of  this  small   (2.25  inches 

O.D.)  toroid  when  full  wound  is  shown  in  Figure  5.1.     The  curve  above 

100  Hz  is  the  manufacturers  specification  which  was  verified  for  a 

typical   full  wound  core  using  a  GR  1650A  inductance  bridge  down  to 

300  Hz.     The  curve  below  100  Hz  is  the  extrapolation  of  the  author. 

The  Q  at  59  Hz  was  checked  using  the  circuit  shown  in  Figure  5.2. 

R    was  determined  by  measuring  the  dissipation  factor,  D,  of  the 

870  uf  capacitor  bank  and  using  the  relation 

Rc  ■  ar  •  <5-6> 

D  was  measured  at  60  Hz  using  a  GR  1611 -A  capacitor  bridge  and  found 
to  be  .0045.  Five  copper  lug  connections  are  included  in  the  tank 
circuit  along  with  thirteen  feet  of  No.  10  AWG  solid  copper  magnet 
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wire  for  a  total  measured  resistance  of  15  mft  using  a  GR  1650-A 

resistance  bridge.  Tne  voltage  waveform  across  the  inductors  was 

observed  with  a  Tektronix  549  storage  oscilloscope.  The  8.4  mh 

inductor  shown  in  the  figure  is  four  2.1  mh  inductors  wound  on  the 

t.l.m.  cores.  The  tank  circuit  was  "pumped  up"  using  the  oscillator 

at  59  Hz  after  which  the  switch  was  opened  and  the  decay  waveform 

observed  on  the  storage  oscilloscope.  Twelve  separate  measurements 

were  made  covering  the  range  of  initial  voltages  from  0.02  v  peak 

to  0.4  v  peak  -  corresponding  to  peak  initial  currents  in  the  range 

from  7  ma  to  129  ma.  In  each  case  the  time  was  measured  for  the 

voltage  across  the  tank  to  fall  to  1^  of  its  initial  value.  The 

times  measured  ranged  from  a  low  of  141  m  sec  at  0.4  v  peak  to 

157  m  sec  below  0.1  v  volt  peak.  The  total  resistance  in  the  tank 

was  determined  from 

R     -  2L       _  16  8  rnh  ,5  7) 

KT  "  time  "  TTrne       *  v°  '' 


idnere 


RT  "  Rc  +  RH  +  Rdc  +  Rac  •  (5"8) 


Finally,  Q  of  the  inductors  was  determined  from 

Q  =  R     o       •  (5-9) 

Rdc   Rac 
The  twelve  Q  values  at  59  Hz  ranged  from  a  low  of  34.6  at  .4  v  peak 

(initial   I  =  129  ma),  to  a  high  of  40  for  initial   voltage  values 

below  0.1   v  peak  (32  ma).     A  similar  "pumping"  technique  was  used 

to  measure  the  inter-winding  capacitance  of  a  single  t.l.m.  inductor 

of  2.1  mh.     A  series  circuit  was  formed  of  an  H.P.  4204A  variable 

frequency  oscillator,  a  4.7  meg  ft  resistor,  and  a  single  2.1  mh 
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t.l.m.   inductor.     The  voltage  across  the  coil  was  observed  while  the 

frequency  was  increased.     The  voltage  peaked  sharply  at  300.0  KHz 

implying  an  inter-winding  capacitance  of  134  p.f. 

The  maximum  inductor  Q  shown  in  Figure  5.1   occurs  when  the 

losses  in  the  core  of  the  inductor  equal   the  losses  in  the  resistance 

of  the  windings.     At  low  frequencies  the  winding  losses  predominate; 

at  maximum  Q  the  two  losses  are  equal,  and  for  frequencies  above 

Q  the  core  losses  predominate.     The  initial   slope  of  the  inductor 

max . 

Q  curve  is  determined  by  winding  losses  which  remain  essentially 

constant  in  this  frequency  range.         In  other  words  -  at  low  frequencies, 

the  Q  is  determined  by  winding  resistance  alone.     In  fact,  a  t.l.m. 

inductor  of  2.1  mh  at  60  Hz  has  wL  =  800  mQ     and  a  winding  resistance 

of  20-mft  indicating  a  Q  of  40.     To  bring  the  Q  of  the  t.l.m.  section 

down  to  20,  for  example,  20  mft  of  resistance  is  to  be  inserted  into 

the  inductor  winding  -  sufficient  to  bring  the  overall  Q  of  the 

section  at  60  Hz  down  to  about  20.     Other,  larger  values  of  resistance 

would  have  to  be  added  to  produce  Q's  less  than  20  for  modeling  lower 

Q  lines.     The  second  curve  in  Figure  5.1   represents  a  typical   Q 

characteristic  of  a  section.     Note  that  the  section  characteristic 

is  lower,  and  the  maximum  is  shifted  to  the  right  -  reflecting  the 

fact  that  core  losses  must  reach  a  higher  value  before  they  can 

equal   the  winding  losses  thus  producing  the  Qmav     condition.     Note 

max . 

also  that  having  shown  a  case  in  which  the  winding  resistance  was 
doubled,  the  initial  slope  of  the  section  Q  curve  is  about  one-half 
of  that  of  the  inductor  alone.     The  third  curve  in  Figure  5.1   is  that 
of  the  typical   765  kV  line  transferred  from  Figure  4.1;  this  curve 
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represents  the  highest  Q  which  must  be  achieved  by  a  section.     It  is 

seen  that  setting  the  Q  of  the  section  at  60  Hz  to  twenty  results 

in  the  section  Q  falling  below  that  of  the  line  above  1700  Hz.     Thus, 

if  Q  is  a  critical   factor  above  1700  Hz  in  a  particular  test  on  a 

typical   765  kV  line,  it  will   be  necessary  to  accept  higher  Q  at  sixty 

Hz  by  setting  less  than  20  mfi  in  series  with  the  model   inductors. 

5.3  Design  of  Model  Self  Inductors 

The  inductance  of  a  toroid  is  given  by 

4N2  u  A  x  10"9 
L  =  q henries,  (5-10) 

where 

N  =  number  of  turns 

u  =  relative  permeability 

A  =  core  cross  sectional   area  (cm) 

D  =  mean  core  diameter  (cm)   . 
For  the  Arnold  core  A-218259-2  used  in  this  design,  up  =  205,  so  that 
L  =  259  mh  per  1000  turns  on  any  full  wound  core.     A  full  wound  core 
is  defined  as  one  in  which  the  winding  occupies  three-fourths  of  the 
available  window  area  with  a  packing  factor  of  0.7.     For  any  other 
full  wound  coil   on  the  A-218259-2, 


N  =  1000  J^g     ;  (5-11) 

where  "L"  is  the  desired  inductance  in  mh.  The  inductances  shown  in 
Table  4.5  require  the  turns  shown  in  Table  5.1.  The  listing  of  turns 
required  in  the  8-mile  sections  is  shown  in  Table  5.2.  To  achieve 
the  Q  characteristic  shown  in  Figure  5.1,  it  is  necessary  to  have  the 
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Table  5.1 
Turns  Required  in  4-Mile  t.l.m.    Inductors 


!!» 

\ 

!!c 

93 

91 

90 

94 

92 

91 

95 

93 

92 

96 

94 

93 

97 

95 

94 

98 

96 
Table  5.2 

95 

Turns 

Required 

in  8-Mile  t.l .m. 

Inductors 

\ 

\ 

!k 

131 

129 

127 

133 

130 

129 

134 

131 

130 

136 

133 

131 

137 

134 

133 

139 

136 

134 

core  full  wound.     In  the  case  of  the  4-mile  inductors  this  works  out 
to  be  No.   10  AWG,  and  for  the  8-mile  sections  -  No.   11  AWG.     For  the 
4-mile  sections,  a  nominal   twenty  feet  of  No.   10  is  required  at  1.0 
mft/ft  for  20  mft  per  inductor.     For  the  8  mile  sections,  a  nominal 
twenty-eight  feet  of  No.   11   is  required  at  1.26  mft/ft  for  35    mft 
per  8-mile  inductor.     In  Tables  5.1   and  5.2,  the  first  five  numbers 
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in  each  column  are  the  turn  numbers  at  which  a  tap  must  be  provided. 

Referring  to  equation  (5-10)  it  is  seen  that  the  inductance 

of  a  toroid  is  independent  of  current  if  u  is  independent  of 

current.  A  graph  of  y  vs.  a-c  flux  density  taken  from  reference  18 

is  reproduced  in  Figure  5.3.  It  is  seen  that  u  remains  within  five 

percent  of  nominal  up  to  4000  gauss.  From  basic  inductor  theory: 

.    X   N^   NBA   NuHA  /,  ,«x 

L  "  I  "  I  "  I  "   I  K°~l£' 

H-dJl  =  NI  (5-13) 

ttDH  =  .MI  (5-14) 

H  =  NI/ttD  (5-15) 


Substitute  (5-15)  into  (5-12): 

NBA       N2yA 
I     "     irD 


B       Ny 
I   "  ttD 


(5-16) 
(5-17) 


I   =  Jjfi  (5-18) 

In  equation  (5-18) 

I  =  complex  amplitude  of  current  (amperes) 
N  =  no.  of  turns  in  winding 

D  =  mean  diameter  of  core  in  meters  =   .0463  m 

2 
B  =  peak  flux  density   (wb/m  ) 

u  =  normal   permeability  =  u  u   . 

For  the  t.l.m.   inductors  at  4000  gauss: 

B  =   .4  wb/m2 

u  =  205(4ttx10'7)^  =  2.58xl0"4 

D  =  0.0463  m 

N  =  93  or  133    . 
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So  that  I  peak  =  T337-  amperes 

=  1 . 70  amperes , 
and  the  inductance  of  the  t.l.m.   inductors  remains  within  five  percent 
of  nominal   up  to  1.70  amperes  r.m.s.,  or  31   p.u. 
5.4     Design  of  Model  Mutual    Inductors 

Figure  2.2  shows  the  mutual   inductors  Lj2(L  J  and  L^o (Lbc) 
as  transformers  between  the  phases.     Since  the  mutual   inductances  are 
small,  and  since  there  is  sufficient  space  remaining  in  the  windows 
of  the  self  inductors,  L  .    is  wound  over  L   ,  and  L.      is  wound  over  L. 
The  mutual   inductances  to  be  provided  in  the  4-mile  sections  are 
shown  in  Table  4.6;  those  for  the  8-mile  sections  are  simply  double 
the  values  in  Table  4.6.     From  basic  mutual   inductor  theory  (See  e.g. 
reference  15  p.p.   178  ff)  it  is  known  that,  for  example, 

Lab  '  kv/tLab>s  La  <5-19> 

Lab  =  tljJ^-  ,  (5-20) 

and  from  equation  (5-11)  for  the  model   core, 


<Nab>s*1000^ir  •  (5-21) 

Where  (all  in  mh): 

L  .  =  mutual  inductance  between  phases  "a"  and  "b" 

(L  .  )  =  self  inductance  of  L  . 

k  =  coefficient  of  coupling  between  L  and  L  . 

L"  =  measured  mutual  inductance  when  L  and  L  . 

a     ab 

are  connected  to  reinforce. 
L"  =  measured  mutual  inductance  when  La  and  L  . 
are  connected  to  oppose. 
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(Nab)s  =  No.   of  turns  to  obtain   (L^L,   (or  L  .). 
Several   tests  were  performed  to  determine  k  by  winding  six,  seven,  and 
eight  turns  of  (L  .  )     over  an  L     inductor  of  2.49  mh.     Knowing  the 
number  of  turns,   (L  .  )     was  obtained  from  (5-21),  and  L  .    was  obtained 
as  indicated  by   (5-20).      (Lat,)s  and  Lah  were  tnen  substituted  into 
(5-19)  to  determine  k.     Under  these  conditions,  k  for  the  t.l.m. 
inductors  was  found  to  be  0.986.     Solving  (5-19)  for  (L  ,)     and 
substituting  into  (5-21), 

<Nab>s*63Lab    17    •  <5-21> 

a 

With  L  taken  as  2.49  mh, 
a 

(Nab>s  *  39-8  Lab>  <5"22> 

and  the  number  of  turns  to  achieve  the  mutual  inductance  values,  L  .  , 

of  Table  4.6  are  6,  7,  8,  9,  10,  and  13,  so  that  taps  must  be  provided 

at  6,  7,  8,  9,  and  10  turns.  With  the  taps  set,  the  self  inductance 

corresponding  to  each  tap  can  be  found  using  equation  (5-21)  and  then 

the  values  of  L  .  with  other  values  of  L  can  be  found  using  (5-19); 

the  results  are  shown  in  Table  5.3.  A  similar  procedure  applied  to 

Lbc  yields, 

<Nbc>s  =  40'5  Lbc 
and  still   the  taps  are  at  6,  7,  8,  9,  and  10  turns,  with  the  13th 

turn  providing  L      -  0.32  mh.     A  list  of  L      values  for  a  4-mile 

section  is  given  in  Table  5.4.     In  order  to  wind  L  .    and  L      over 

L    and  L    respectively,  it  is  necessary  to  use  four  feet  of  No.   12 

AWG  magnet  wire  at  1.59  mft/ft  for  an  additional   resistance  of  6.4  mQ 

in  series  with  the  phase  "b"  and  "c"  lines.     A  list  of  mutual 

inductance  values  to  be  built  into  the  8-mile  sections  is  shown  in 
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Table  5.5. 


Table  5.5 

Lab(mh) 

Lbc(mh) 

0.30 

0.30 

0.36 

0.36 

0.40 

0.40 

0.46 

0.46 

0.50 

0.50 

0.66 

0.56 

Table  of  Mutual    Inductance  Values 
Designed  into  the  8-Mile  Sections 

The  taps  to  yield  (Lab)8  (using  La  =  4.98)  are  8,  10,  11,  13,  14,  and 

18;  taps  for  (Lbc)8  (using  Lfa  =  4.82)  are  8,  10,  11,  13,  14,  and  17. 

The  corresponding  values  of  (L  .  )     for  the  8-mile  sections  are  (in  mh): 

0.0166,  0.0259,  0.0312,  0.0435,  0.0505,  and  0.0834.     For  (Lbc)s: 

0.0166,  0.0259,  0.0312,  0.0435,  0.0505,  and  0.0745.     These  tabulations, 

along  with  those  of  Tables  5.3  and  5.4  show  that  the  self  inductance 

of  a  mutual  winding  is  negligible  when  compared  to  the  self  inductance 

of  the  line.     In  addition,  Tables  5.3  and  5.4  indicate  that  the  change 

in  L  .    or  L      resulting  from  a  change  in  L     or  L    over  their  entire 

range  is  less  than  five  percent  and  therefore  negligible  also.     In 

order  to  wind  L  .    and  L,  „  over  L,  and  L    in  the  8-mile  sections,  it 
ab  be  a  b 

is  necessary  to  use  approximately  six  feet  of  No.   12  AWG  magnet  wire 
at  1.59  mft/ft  for  an  additional   resistance  of  9.5  mft  in  series  with 
the  phase  "b"  and  "c"   lines. 
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CHAPTER  6 
CAPACITOR  SELECTION 
6.1  Factors  Determining  Capacitor  Selection 

As  discussed  in  detail  in  Chapter  4,  eight  diverse  transmission 
line  types  were  analyzed  in  the  EPSEL  Line  Parameters  Program,  and 
these  eight  types  were  taken  to  cover  the  range  of  Typical  EHV  lines. 
The  model  capacitance  values  implied  by  these  lines  for  a  4-mile  model 
section  are  shown  in  Tables  4.4a  and  4.4b.  A  complete  listing  of  the 
capacitance  parameters  to  be  covered  by  both  the  4  and  8-mile  model 
sections  is  shown  in  Table  6.1.  Having  scaled  the  eight  representative 
lines  to  result  in  a  minimum  inductance  range,  the  range  of  each 
capacitance  is  seen  to  be  somewhat  larger.  It  is  proposed  to  cover 
the  range  of  each  parameter  of  Table  6.1  with  a  bank  of  six  capacitors 
wired  and  brought  to  the  operating  panel  so  that  any  combination  of  the 
six  capacitors  in  each  bank  can  be  arranged  in  parallel.  Thus,  with 
one  or  two  capacitors  to  set  the  lowest  value  in  a  range,  approximately 
five  or  four  capacitors  are  available  to  cover  the  remainder  of  the 
range.  Thus  there  are  a  maximum  of  120  and  minimum  of  24  capacitance 
values  available  to  cover  the  range  of  each  capacitance  parameter 
indicated  in  Table  6.1.  With  this  "dense"  provision  of  model  values 
it  is  possible  to  set  any  value  in  the  span  within  a  very  few  percent 
(^3%)  if  the  tolerance  of  the  capacitors  is  less  than  a  few  percent. 
It  is  proposed  to  provide  a  two  to  three  percent  tolerance  on  any 
capacitance  selected  by  the  operator  so  that  the  over-all  capacitance 
provision  of  the  model  is  within  about  5%  of  the  value  required  to 
model  any  given  line.  With  expected  transient  overvoltage  maxima 
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of  about  four  p.u.,  all  capacitors  must  be  capable  of  withstanding  short 
duration  peaks  up  to  500  or  600  volts.  In  addition  to  tolerance  and 
voltage  requirements  the  capacitors  must  have  extremely  low  losses. 

An  indication  of  the  dissipation  factor  required  of  the  model 
capacitors  can  be  obtained  by  a  simple  analysis  of  the  "dissipation 
factor"  of  a  typical  line.  Considering  the  positive  sequence  parameters 
of  a  type  41 A  line  using  the  values  given  in  Table  4.1, 

BP  S 
CD  c  =  -^^-  =  .0181  uf/mi. 

Cp  s  =  .0725  yf/4  mi. 

Rp  s  =  .216  n/4  mi. 
Using  the  "series  model"  for  a  lossy  capacitor: 

0pS>/4  mi  =  wRptS#Cp>s<  (6-1) 

0p  s  =  0.0000059/4  mi.  Q   60  Hz  (6-2) 

Op  s     =  0.0000985/4  mi.       <a  1000  Hz  (6-3) 

That  the  values  given  by  equations  (6-2)  and  (6-3)  are  the  values 
required  of  a  4-mile  model  section  can  be  seen  by  a  corresponding 
calculation  for  the  model.     Using  Table  4.4a  for  the  type  41A  model: 

CP.S.  ■  CJ  +  3CM  ■  Cd"  *  V  -  CC"  +  Cab  +  Cac  +  Sc 

Cp  s     =   .230  uf/4  mi. 

T       R     +  R.    +  R 
Rp  s     =  R^  =  — ^ 2-  =   .068  fi/4mi. 

Up  s     =  0.0000059/4  mi.        @  60  Hz  (6-4) 

Up  s     =  0.0000985/4  mi.        0  1000  Hz.  (6-5) 
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at  low  currents.  The  basic  relation  governing  the  Q  of  an  inductor  can 
be  written  as  (for  the  frequencies  of  interest  here) 

Q  =  — 4^ o-  t5-1* 

r  .  r  +  r  r 

dc     ac 
Where 

L  =  Inductance  in  henries 

co  =  2tt  x  frequency 

R,     =  dc  resistance  of  the  winding   (ft) 

Ra„  =  equivalent  core  loss  resistance  (ft) 

2 

With  these  definitions  it  is  seen  that  R  I  corresponds  to  both  the 

ac  - — ■— 

hysteresis  and  eddy  current  losses  (See  reference  17  p.p.  58ff.)  and 

thus  represents  the  total  core  loss.  R_.  the  core  loss  coefficient, 

ac 

18 
is  itself  a  function  of  current.     For  powdered  cores  the  manufacturer 

provides  data  on  R      explicitly,  whereas  for  Fe-Si   the  entire  term 

ac 

R     I     is  provided.     If  the  data  provided  by  the  manufacturer  is 
ac 

converted  to  the  terms  of  equation  (5-1),  it  is  found  that  for  Fe-Si 
the  core  loss,  W,   is  given  by  (approximately) 

W  =  k^uOl1-7,  (5-2) 

whereas  for  the  powdered  Mo-Ni   the  loss  is  given  by   (approximately) 

W  =  RacI2  =  k2(w)I2*4  .  (5-3) 

In  equations   (5-2)  and  (5-3),  k,(to)  and  ko((jo)  are  core  loss  constants 
and  are  a  function  of  frequency.     If  equations   (5-2)  and  (5-3)  are 
substituted  in  equation  (5-1)  there  results,  respectively, 

0  = wkl (5.4) 

gFe-Si       I     r2   .   .    ,   wl.7     •  ^  q) 

RdcI     +  k^wjl 

and 
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occurs  primarily  at  low  frequencies  in  connection  with  switching  opera- 
tions which  connect  the  terminals  of  a  charged  line  to  large  inductive 
reactances.  Natural  frequencies  of  the  order  of  40  Hz  are  known  to  occur 
under  these  conditions  with  over-all  Q  approaching  300. 

Most  capacitor  manufacturers  specify  or  test  dissipation  factor  at 
1000  Hz,  so  that  at  1000  Hz  it  is  possible  to  measure  capacitance  and 
compute  resistance  using  either  the  parallel  or  series  model  of  a  lossy 
capacitor: 

D  =  uRscs  ■  srJc^  ■  <6-6> 

The  values  of  R  or  R  so  obtained,  however,  are  true  only  for  a  limited 
range  of  frequencies  and  cannot  be  used  at  40  Hz.  As  frequency  decreases 
toward  zero  the  value  of  R  (resistance  in  the  parallel  model)  approaches 

the  insulation  resistance  of  the  dielectric  -  which  for  polystyrene  is 

20 
of  the  order  of  ten  million  megohms  for  a  0.15  uf  capacitor.    The 

dissipation  factor,  D,  of  a  typical  .15  yf  polystyrene  capacitor  at 

20 
1000  Hz  is  given  as  0.0001.    This  value  of  D  implies  R  =  10.6  megohms 

at  1000  Hz.  Few  materials  other  than  Teflon  or  a  modified  polystyrene 

can  equal  this  low  value  of  D  or  high  value  of  R  at  1000  Hz.  If 

equation  (6-6)  is  solved  for  R  in  terms  of  w,  R$,  C  ,  and  C 

(C  =  C  ),  and  R  set  equal  to  1.0  mft,  C  =  0.15  uf,  and  f  =  40  Hz  the 
s    p       s  s 

value  of  R  obtained  is  4.4  million  megohms.  Thus  in  order  to  obtain 
R  as  low  as  .001  ft,  it  is  necessary  to  have  R  as  high  as  4.4  million 

megohms  at  40  Hz.  Little  data  exists  on  any  dielectric  material  below 

21 

100  Hz  -  standard  reference  works   supply  no  data  below  100  Hz  on  any 

material  and  give  no  information  vs.  frequency  on  polystyrene.  In  an 
effort  to  obtain  an  estimate  of  R  for  the  0.15  uf  polystyrene  capacitor 
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at  40  Hz,  it  was  assumed  that  R  remains  constant  at  the  10.6  megohm 
value  down  to  100  Hz  and  then  rises  linearly  to  10  million  megohms  at 
d.c.  This  procedure  gives  R  as  6  million  megohms  at  40  Hz.  (This  is 
a  conservative  estimate  in  the  opinion  of  the  author.)  No  material 
less  costly  than  polystyrene  can  approach  within  an  order  of  magnitude 
of  this  value.  Polystyrene  capacitors  are  therefore  selected  to  form 
the  "base"  capacitance  of  each  capacitor  bank  in  the  model. 
6.2  Capacitors  Selected  for  Model  Sections 

Since  the  polystyrene  capacitors  conservativety  (but  not  by 
orders  of  magnitude)  meet  all  line  requirements,  a  compromise  in  favor 
of  lower  cost  is  made  in  selecting  the  trimming  capacitors  in  each 
bank.  Thus,  while  polystyrene  capacitors  will  be  made  to  provide 
80  percent  (approximately)  of  any  capacitor  value  selectable  by  the 
operator,  the  remaining  20  percent  (approximately)  is  to  be  made  up  of 
lower  cost  mylar  or  silvered  mica  capacitors.  The  latter  capacitors 
have  losses  ten  to  fifty  times  higher  than  the  polystyrene  but  they 
constitute  only  ten  to  twenty  percent  of  a  selected  capacitor  value 
and  so  should  not  appreciably  lower  the  high  quality  provided  by  the 
polystyrene  base  capacitors.  The  severe  low-loss  constraint  applies 
mainly  to  the  765  kV  line,  so  that  the  polystyrene  capacitors  in  each 
bank  are  selected  so  that  for  the  765  kV  case  they  will  constitute  no 
less  than  90%  of  the  required  capacitance  of  a  4-mile  section.  The 
polystyrene  capacitors  selected  for  the  t.l.m.  are  Arco  Manufacturing 
Co.,  type  PJ,  400  VDC,  +  2%,  tubular  capacitors.  The  400  V  units 
were  judged  adequate  since  each  capacitor  receives  a  minimum  of  two 

sixteen  hour  burn-in  periods  with  300  percent  of  rated  voltage  applied 

20 
at  85°C.    The  mylar  capacitors  are  Cornell  Dubilier  Co.,  type  WMF6, 
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600  VDC,  +  10  percent,  tubular  capacitors.  Two  small  capacitance  values 
below  the  range  of  standard  mylar  sizes  (330  pf  and  680  pf)  are  Cornell 
Dubilier  Co.,  type  CD19FD,  silvered  mica,  500  VDC,  +  5%.  The  capacitor 
types  and  sizes  selected  to  span  each  capacitance  parameter  (the 
capacitor  banks)  are  shown  in  Table  6.2. 
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CHAPTER  7 
RESISTANCE  AND  LOSS  CONSIDERATIONS 
7.1  Over-all  Q  Considerations 

As  discussed  in  Chapters  4  and  5,  the  Q  requirements  of  this 
t.l.m.   are  dictated  by  the  need  to  model   the  765  kV  line.     This  con- 
straint is  particularly  severe  in  the  4-mile  model   sections  since 
connection  resistance  is  expected  to  take  up  any  margin  provided  by 
designing  the  inductors  with  double  the  required  Q  at  60  Hz.     In 
Chapter  4  it  is  indicated  that  the  positive  sequence  Q  of  a  typical 
765  kV  line  is  twenty.     The  resistance  required  in  each  phase  to  achieve 
this  result  is  shown  in  Table  4.4a  -  37  mft  in  phases   "a"  and  "c"  and 
43  mft  in  phase  "b".     As  discussed  in  Chapter  5,  20  mft  is   "provided"  in 
each  phase  by  the  resistance  of  the  wire  to  wind  the  self  inductors 
L   ,  L,  ,  and  L  .     An  additional   6  mfl  is  provided  in  phases   "b"  and  "c" 
due  to  the  fact  that  these  phases  contain  the  resistance  of  the  windings 
for  L  .    and  L      respectively.      (Actually,  L  .    and  L       for  the  typical 
765  kV  line  are  provided  by  the  7- turn- tap  so  that  the  resistance 
added  by  the  L  .    and  L      windings  will   be  approximately  4  mft  in  the 
4-mile  sections.     The  operating  panel   resistance  to  be  provided  at  this 
point  in  the  accounting  is  17  mi^  in  phase  "a",  19  mft  in  phase  "b",  and 
13  mft  in  phase  "c".     But  connection  resistance  has  not  yet  been 
accounted  for  in  these  figures. 

In  a  typical  model  set  up  (See  Chapter  8),  six  to  eight  banana 
jack  connections  will  be  made  in  series  in  each  section  from  input  to 
output,  so  that  the  resistance  of  these  must  be  accounted  for.  The 
jacks  selected  for  this  design  are  E.F.  Johnson  Co.,  standard  banana 
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jack,  cadmium  plated  insert,  catalog  no.   108-0900-001.     The  corresponding 
plug  is  E.F.  Johnson  Co.,  standard  plug,  catalog  no.   108-0300-001.     Six 
of  these  jacks  and  plugs  were  connected  in  series  as  shown  in  Figure  7.1. 
Sixteen  inches  of  solid,  no.   14  AWG,  copper  bus  bar  wire  was  used,  and 
all   connections  were  heavily  soldered.     The  total   resistance  was  measured 
at  12  mft  with  the  GR  1650A  resistance  bridge.     The  resistance  of  the 
no.   14  wire  is  3.3  mft  so  that  the  resistance  of  each  connection  averages 
to  about  1.5  mft.     All   plugs  and  jacks  were  new,  clean,  tight  fitting, 
and  all   connections  were  made  and  remade  several   times  -  always  with  the 
same  result.     The  resistance  accounting  to  this  point  is: 
R     =  8  mft,   R,    =  5  mft,   R     =  1  mft.     In  other  words,  at  this  point  the 
resistance  values  required  at  the  front  panel   to  bring  the  section 
positive  sequence  Q  down  to  twenty  at  60  Hz  are  the  values  cited  just 
above.     Each  4-mile  section;  exclusive  of  external   hookup  wire,  but 
including  operator  panel   connector  resistance  is  seen  to  have  built 
into  it, 

Ra(antrinsic)4  =  29  mft  (7-1) 

Rb(intrinsic)4  =  38  mft  (7-2) 

R  (intrinsic),  =  38  mft  (7-3) 

To  provide  an  exact  positive  sequence  Q  of  twenty  at  60  Hz,  it  is 
necessary  to  provide  for  operator  selection;   R     =  8  mft,  R,    =  5  mft, 
and  R     =  1  mft.     As  pointed  out  in  Chapter  5  in  connection  with 
Figure  5.1   however;  setting  the  Q  to  twenty  at  60  Hz  for  a  765  kV  line 
model  will   result  in  the  model   Q  falling  below  line  Q  above  about 
1700  Hz.     Thus,  it  may  be  desirable  under  certain  circumstances  to  set 
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Figure  7.1 
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Q  at  60  Hz  slightly  above  twenty  to  ensure  that  the  model   Q  remains 

above  line  Q  up  to  3060  Hz. 

An  accounting  procedure  identical   to  that  above  applied  to 

the  8-mile  model  sections  gives  the  results, 

Ra(intrinsic)8  =  44  mft  (7-4) 

R,  (intrinsic)o  =  57  ma  (7-5) 

R  (intrinsic)g  =  57  ma  (7-6) 

The  values  of  R3,  R,  ,  and  R     required  in  the  8-mile  sections  to  model   a 
a       d  c 

765  kV  line  are   (doubling  the  appropriate  values  in  Table  4.4a): 

(Ra)765  -  74  ma  (7-7) 

(Rb)765  -  86  m®  (7-8) 

(Rc)765  =  74  ma  (7-9) 

It  is  seen  that  the  8-mile  sections  have  considerably  more  margin  than 
the  4-mile  sections  and,  in  fact,  using  the  intrinsic  values  of 
equations   (7-4)- (7-6),   the  positive  sequence  Q  of  the  8-mile  sections 
can  be  set  as  high  as  31   at  60  Hz. 
7.2     Resistance  Selection 

Referring  to  Table  4.2   (excluding  the  1500  kV  lines)   the  values 
of  positive  sequence  Q  to  be  provided  are  seen  to  range  from  about 
eight  to  twenty.     To  achieve  the  Q  values  lower  than  twenty,  resistance 
values  must  be  provided  at  the  operating  panel    for  insertion  with  the 
intrinsic  resistance  of  the  section.     For  the  4-mile  sections,  the 
lowest  resistance  value  to  be  made  available  is  clearly  zero  since 
this  is  the  required  value  to  ensure  positive-sequence  Q's  just  over 
twenty.     For  tiie  8-mile  sections  a  similar  zero  base  will   be  provided 
since  it  would  simply  not  make  sense  to  deny  the  operator  the  selection 
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of  positive  sequences  Q's  to  thirty  should  the  occasion  arise.     A  list- 
ing of  Q  values  to  be  provided  for  each  phase  is  shown  in  Table  7-1. 

Table  7.1 
4-Mile  Section  8-Mile  Section 


% 

% 

^c 

^a 

% 

% 

30 

22 

21 

40 

30 

28 

25 

19 

18 

34 

25 

24 

20 

16 

15 

28 

20 

20 

15 

13 

12 

22 

16 

16 

12 

11 

10 

16 

12 

11 

10 

9 

8 

11 

9 

8 

Table  of  Q  Values  Provided 
in  the  T.L.M.  Sections 

The  operating  panel   resistance  values  to  achieve  the  Q's  shown  in 
Table  7.1   are  shown  in  Table  7.2. 

Table  7.2 
4-Mile  Sections  8-Mile  Sections 


Ra 

Rb 

Rc 

R 
a 

Rb 

Rc 

0 

0 

0 

0 

0 

0 

6 

7 

8 

8 

11 

12 

15 

15 

17 

19 

28 

25 

30 

27 

30 

36 

49 

46 

44 

40 

43 

66 

85 

92 

59 

67 

64 

116 

132 

147 

Table  of  Operating  Panel   Resistance  Values 
Corresponding  to  Q's  of  Table  7.1      (R  in  mfl) 
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The  resistance  values  indicated  in  Table  7.2  were  obtained  using 

w(L.)d  o 
R.  = q-1^  -  R.(intrinsic)  (7-10) 

J 

(La)4  =  2.34  mh,  (Lb)4  =  2.26  mh,  (Lc)4  =  2.18  mh 
(La)8  =  4.68  mh,  (Lb)8  =  4.52  mh,  (Lc)8  =  4.36  mh 
i  =  a,b,c  for  each  section  length 
Q.  is  the  corresponding  Q  value  from  Table  7.1. 

<J 

It  is  seen  in  practice,  rather  than  setting  positive  sequence  Q  per  se, 
the  Q  of  each  phase  is  set  individually  using  Tables  4.4,  7.1,  and  7.2. 
The  resistance  indicated  in  Table  7.2  is  achieved  by  connecting 
resistances  in  series  so  that  each  indicated  value  is  the  sum  of  the 
resistances  required  to  form  each  of  the  preceeding  resistors.  The 
actual  resistors  required  to  obtain  the  operating  panel  values  indicated 
in  Table  7.2  are  shown  in  Table  7.3. 

Table  7.3 

4-Mile  Sections  8-Mile  Sections 

R   *s       R  ,'s       R   's       R   's       R  ,'s       R   's 
a  b  c  a  b  c 

6  7  8  8  11  12 

9  8  9  11  17  13 

15  12  13  17  21  21 

14  13  13  30  36  46 

15  27  21  50  47  55 

Table  of  Resistors  in  the  T.L.M.  Sections 

The  resistors   indicated  in  Table  7.3  are  to  be  "manufactured"  by 
cutting  appropriate  lengths  of  nickel   or  copper-nickel   alloy   ("cupron") 
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resistance  wire.     A  very  large  selection  of  resistance  wire  sizes  and 
types  are  available  from  the  Driver  Harris  Co.   of  .Newark,  N.J.,  a  type 
particularly  suited  to  the  smaller  resistors  above  is  the  company's 
0.035"  diameter  nickel  wire  of  0.045  ohms  per  foot.     For  the  larger 
resistors  no.   20  AWG  "Cupron"   (copper-nickel)  wire  with  a  nominal 
resistance  of  0.287  ohms  per  foot  can  be  used.     The  latter  wire  can 
safely  carry  a  direct  current  of  three  amperes   indefinitely  but  is 
hot  to  touch  at  this  current  level.     Three  rms  amperes   is  far  in 
excess  of  that  required  to  simulate  the  most  severe  faults,  so  that 
no.   20  AWG  resistance  wire  should  be  adequate  for  this  t.l.m. 
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CHAPTER  8 
MODEL  PHYSICAL  CONSTRUCTION 

8.1  Electrical   Construction  of  a  Typical  Section 

As  discussed  in  Chapter  4,  this  t.l.m.   consists  of  thirty  4-mile 
sections  and  thirty-five  eight  mile  sections.     A  simplified  schematic 
diagram  of  the  wiring  of  a  typical   section  is  shown  in  Figure  8.1   and 
is  the  same  for  both  the  4-mile  and  8-mile  sections.     The  component 
values  are  not  indicated  since  they  are  different  for  the  two  section 
types,  but  the  required  values  can  be  taken  from  the  tables  of  Chapters 
5,  6,  and  7.     In  wiring  each  section,  all   hookup  wire  in  a  direct  path 
from  input  to  output  is  No.   10-stranded,  tinned,  copper  wire.     Capacitor 
hookup  is  effected  with  No.   16  stranded  copper  wire.     All   connections 
are  soldered. 

8.2  T.L.M.   Physical   Construction 

Each  t.l.m.  section  is  constructed  on  the  chassis  shown  in 
Figure  8.2.     The  sixty-five  chassis  are  made  of  plexiglass  of  the  sizes 
indicated  -  the  operator's  panel   is  of  black  plexiglass.     The  layout 
of  the       operator's  panel   is  shown  in  Figure  8.3.     The  jacks  are  color 
coded  by  "phase"  and  the  circuit  diagram  is  painted  by  a  silkscreen 
stenciling  process.     Each  of  the  sixty  five  chassis   is  assembled  in  a 
larger  console  as  shown  in  Figure  8.4.     The  chassis  slide  into  the 
console  on  rails  and  are  secured  with  screws  through  the  operating 
panel.     The  chassis  are  inserted  into  the  console  from  both  "front" 
and  "rear",  so  that  twenty-four  sections  are  accommodated  per  console. 
Three  consoles  such  as  that  shown  in  Figure  8.4  are  required  with  a 
fourth   "half-console"  for  the  output  section  and  terminal   equipment. 
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CHAPTER  9 
SUMMARY  AND  RECOMMENDATIONS 
9.1  Summary 

This  report  presents  in  detail   the  considerations  in  the  design 
of  a  transmission  line  model   for  use  in  the  M.I.T.   Transmission  System 
Simulator.     The  design  includes  the  modeling  of  the  inductive  effects  of 
non-transposition  by  incorporating  two  linear  mutual   inductors  between 
phases.     The  inter-phase  resistance  effects  are  approximated  by  assuming 
that  the  three  inter-phase  resistances  are  equal.     A  circuit  diagram  of 
one  three-phase  L-section  is  shown  in  Figures  2.2  and  8.1.     Note  that 
the  mutual   inductors  are  physically  wound  on  the  respective  self  inductor 
core  and  that  the  return  circuit     is  provisional   pending  a  complete 
design  due  in  September,  1971. 

In  order  to  deduce  the  number  of  sections  required  for  the  model 
and  establish  a  quantitative  measure  for  predicting  the  performance  of 
model   types  with  various  terminations,  the  pole  error  analysis  of  Chapter 
3  was  developed.     This  analysis  is  based  upon  a  comparison  of  the  poles 
of  the  driving  point  impedance  of  actual   lines  with  those  of  proposed 
models.     Solutions  for  the  poles  of  tt  and  L-section  models  were  obtained 
by  computer  for  models  up  to  forty  sections  with  open  and  short  circuited 
terminations.     The  pole  error  curves  of  Chapter  3  show  in  quantitative 
terms  the  convergence  of  the  poles  of  various  models   toward  those  of 
continuous  lines  with  similar  terminations.     A  technique  is  outlined 
(See  Figure  3.11)  in  which  the  pole  error  curves  are  used  to  deduce 
the  number  of  sections  required  to  meet  any  arbitrary  accuracy  spec- 
ification.    The  specification  set  for  the  M.I.T.  model   is  that  the 
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error  in  the  pole  frequency  of  the  model  shall  be  within  five  percent 
of  the  frequency  of  the  corresponding  actual  pole  over  the  range  from 
sixty  to  3060  Hz  for  line  lengths  of  sixty  to  four  hundred  miles  for 
both  open  and  short  circuited  terminations.  This  criteria  led  to  the 
decision  to  construct  thirty  four-mile  model  sections  and  thirty- five 
eight-mile  sections   ("L"  configuration). 

In  Chapter  4  parameter  data  for  eight  representative  EHV  line 
configuration  types  are  presented  in  order  to  deduce  the  parameter  values 
required  in  the  four  and  eight-mile  model  sections.     Several   scaling 
techniques  are  suggested  for  scaling  the  parameters  of  the  diverse  line 
types  to  the  parameters  provided  in  the  model.     Careful   scaling  suggests 
that  only  a     narrow  range  of  component  values  need  be  provided  for  each 
parameter  in  the  model.     In  the  M.I.T.  design,  six  values  are  provided 
for  each  of  the  sixteen  parameters  to  be  modeled. 

The  quality  factor  (Q)  of  typical   EHV  lines  are  seen   (Figure  4.1) 
to  vary  markedly  with  frequency  in  the  range  from  sixty  to  3060  Hz. 
Positive  sequence  Q  increases  monotonically  to  over  175  for  the  highest 
Q  lines  studied,  and  zero-sequence  Q  maximizes  near  1000  Hz  and 
generally  remains  constant  up  to  3060  Hz.     The  requirement  to  model 
these  Q  characteristics  is  a  major  factor  in  the  design  of  all  model 
inductors  and  in  the  selection  of  all  model   capacitors.     The  positive 
sequence  Q  characteristics  is  achieved  by  the  use  of  molybdenum- 
nickel   toroids  for  all   phase  inductors.     In  fact,  the  model   inductors 
have  an  intrinsic  Q  characteristic  approximately  double  the  value 
required  to  model   the  highest  Q  EHV  lines  in  operation.     Preliminary 
plans  for  modeling  the  return  path  or  zero-sequence  inductance  and 
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Q  suggest  a  separate  modeling  of  the  ground  wire  path  and  earth  return 
path  -  the  final   design  is  to  be  completed  by  September,  1971. 

A  careful   account  of  wire  and  connector  resistance  shows  that 
for  the  section  lengths  chosen,  every  mi  Hi  ohm  significantly  reduces 
the  intrinsic  Q  provided  by  the  model   inductors.     To  avoid  excessive 
losses  and  provide  for  operator  adjustment  of  Q,  all   connections  are 
soldered,  and  all  wire  on  a  direct  path  from  input  to  output  is 
no.   10  AWG  hook  up  wire.     To  avoid  further  losses,  and  to  model   the 
condition  of  a  charged  line  isolated  with  a  high  Q  reactor  at  its 
terminals,  extremely  low  loss  polystyrene  capacitors  are  incorporated 
into  each  t.l.m.   capacitor  bank.     Approximately  eighty-five  percent 
of  any  model   capacitance  parameter  selectable  by  the  operator  is  low 
loss  polystyrene  in  parallel  with  approximately  fifteen  percent 
additional    "trimming"  capacity  provided  by  mylar  units.     Overall 
capacitance  tolerance  in  the  t.l.m.   is  within  five  percent  of  any  value 
selected.     Six  capacitors  are  provided  for  each  capacitance  parameter 
and  all   capacitors  are  available  for  operator  connection  into  any 
desired  configuration. 

The  final   result  for  the  M.I.T.  t.l.m.   is  a  total  of  sixty-five 
L-sections  for  modeling  line  lengths  to  400  miles.     Average  Q  in  the 
four-mile  model   sections  is  a  maximum  of  24  at  60  Hz;  the  maximum 
average  in  the  eight  mile  sections  is  32.     All   capacitance  and 
inductance  selection  within  the  range  of  a  parameter  is  within  five 
percent  of  any  desired  value.     A  similar  five  percent  tolerance  is 
planned  for  the  return  path  inductance  and  Q.     Q  values  less  than  the 
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maxima  given  above  are  obtained  by  the  insertion  of  appropriate  resistors 
provided  in  each  section. 
9.2     Recommendations 

The  design  of  this  t.l.m.   has  progressed  in  parallel  with  its 
actual   construction  and  component  procurement.     This  was  attempted  in 
view  of  a  master's  thesis  time  constraint.     The  basic  scale  factor 
between  the  model   and  lines  was  selected  in  order  to  make  the  t.l.m. 
compatible  with  existing  laboratory  equipment.     In  fact,  the  t.l.m. 
ultimately  will   represent  an  investment  several   times  the  cost  of  the 
TSS  transformer  which  dictated  the  scale  ratio.     This  scale  ratio 
(the  0.355  for  the  765  kV  line)  was  selected  and  found  acceptable 
primarily  on  the  basis   that  it  resulted  in  a  feasible,  low  cost  inductor 
core  without  scaling  the  capacitors  out  of  the  range  of  economical 
values  provided  "off  the  shelf"  by  manufacturers.     Not  until   after  the 
cores  had  been  purchased  was  the  significance  of  capacitor  dissipation 
factor  fully  appreciated.     Although  the  required  capacitor  values  are 
well  within  the  range  of  normal   capacitors,  they  are  among  the  highest 
and  most  costly  polystyrene  units.     The  result  was  the  requirement  to 
purchase  capacitors  nearly  three  times  the  cost  of  inductor  cores. 
Polystyrene  capacitors  are  required  regardless  of  what  reasonable  scale 
ratio  is  selected,  but  it  is  suggested  that  perhaps  a  better  cost  trade 
off  could  have  been  achieved  had  the  TSS  transformers  been  redesigned  to 
permit  a  scale  ratio  closer  to  unity.     This  would  require  larger,  more 
expensive,  cores  for  the  t.l.m.   inductors  but  would  significantly  reduce 
the  size  of  capacitors   required.     The  larger  polystyrene  units  are  very 
expensive,  so  that  even  a  small   reduction  in  values  required  could  result 
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in  considerable  savings.  This  savings  could  be  directed  toward  purchasing 
a  better  (larger)  inductor  core  permitting  even  more  margin  in  Q  by  making 
the  connector  resistance  a  smaller  fraction  of  the  total  resistance  of  the 
section.   It  is  recommended  that,  should  a  similar  project  be  undertaken 
in  the  future,  the  entire  design  should  be  completed  in  every  detail,  a 
prototype  section  constructed,  and  a  cost  analysis  completed  before  a 
single  large  scale  purchase  is  made.  With  a  better  cost  trade  off  between 
t.l.m.  components  and  a  more  complete  view  of  overall  t.l.m.  investment 
compared  to  other  TSS  components,  a  more  informed  decision  can  be  made  as 
to  what  part  of  the  TSS  should  be  scaled  to  achieve  compatibility  between 
them  al 1 . 

Concerning  the  pole  error  analysis  of  Chapter  3,  several  comments 
seem  appropriate.  It  should  be  noted  that  none  of  the  pole  error  analysis 
was  carried  to  more  than  39  sections.  This  is  due  to  the  fact  that 
none  of  the  root  subroutines  available  was  able  to  generate  stable 
solutions  beyond  this  point.  Several  algebraic  root  finding  methods  were 
used  in  an  effort  to  obtain  solutions  above  the  40th  degree,  but  none 
succeeded.  This  is  believed  due  to  the  fact  that  all  the  roots  are  in 
the  region  between  minus  one  and  plus  one,  and  as  the  degree  of  the 
polynomial  gets  higher  the  roots  tend  to  crowd  near  the  extremes.  Scal- 
ing of  the  coefficients  of  the  polynomials  was  attempted  but  this 
technique  resulted  in  instabilities  at  even  lower  degree.  Three  of  the 
four  subroutines  used  in  this  problem  are  available  in  the  IBM  System/360 
Scientific  Subroutine  Package  (360-D1-03X  version  III)  using  three 
different  (but  algebraic)  techniques.  The  subroutine  POLRT  uses  the 
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Newton-Raphson  iterative  technique  and  is  explicitly  limited  to  poly- 
nomials of  degree  thirty-six  or  less.     The  subroutine  DPRQD  uses  the 
quotient-difference-displacement  method  -  this  routine  generally 
became  unstable  at  about  the  30th  degree  for  this  problem.     The  DPRBM 
subroutine  uses  the  Bairstow  iterative  method  of  quadratic  factoriza- 
tion;  this  routine  was  the  most  stable  of  the  four  methods  -  giving 
stable  solutions  through  degree  thirty-nine.     A  fourth  subroutine, 
ROOT  1,  using  Muller's  method   (written  by  the  computation  group  at 
H.I.T. )  became  unstable  at  about  degree  twenty-four.     The  Bairstow 
method  is   the  routine  selected  for  all   the  programs  of  Appendix  D. 
No  extensive  effort  was  made  to  determine  exactly  why  the  subroutines 
became  unstable,  nor  were  more  elegant  matrix  methods  attempted.     It 
was  felt  that  solutions  to  degree  forty  gave  information  sufficient 
for  the  purposes  of  this  design,  so  that  no  further  effort  was  expended 
in  this  direction.     Nevertheless  the  instability  problem  is  at  least 
of  academic  value  and  worth  further  study. 
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APPENDIX  A 
MODEL  TRANSFORMATIONS 

A.I     Equations  of  the  General   Non-Transposed  Model 

In  the  general,  non-transposed  model  of  Figure  2.1   Chapter  2, 

the  differential   equations  relating  the  line  to  ground  voltages  and 

the  currents  are  given  by  Kirchoff's  Voltage  Law   (The  Voltage  Law): 

9V_  91  31, 

-JT-   ("llVa  +   <L11+W   3tT  +  Wb  +  <L12+L00>  3t^  + 

•Wc  +  (Li3+Loo>  Sr  (/Vla) 

3V,  31  31. 

"  ZT  '  Wa  +   (L12+W   3T  +<R22+R00>Ib  +   {L22+W  W~  + 

Wc  +   (L23+L00>  If  (A"lb> 

3Vc  3Ia  3lb 

"   3T  =   R00!a  +   (L13+L00)   3t~  +  R00!b  +   (L23+L00}   3t~  + 

(R33+R00^c  +   (L33+lW  JT  (A"lc) 

Similarly,  the  currents  can  be  related  to  the  line  to  ground  voltages  by 
Kirchoff's  Current  Law   (The  Current  Law): 


31                                        3V                  3V,                  3V 

=   (C     +C     +C     )     ■          c                  c 
3z          ^10  L12  L13;    3t         U12  3t         L13   3t 

(A-2a) 

91,                     3V                                       3V.                 3V„ 

-        C                +    (C     +C     +C      }            -   C 
3z             °12  3t         lt20  L12  L23;   3t         L23  3t 

(A-2b) 

31                      3V                  3V,                                      3Vr 

■  -    -  -  c         a  -  r            +  (c    +c    +c    ) 

3t             U13  3t         ^23  3t          ^L30  L13  L23;   3t 

(A-2c) 

A. 2     Transformation  of  the  General   Non-Transposed  Model 

to  a 

"Ba" 

lanced 

Mutual"  Model 

Equations    (A-l )  describing  Figure  2.1   can  be  transformed  to  a 
set  describing  Figure  2.2  using  simple  algebra.     Let  Ia+Ik+Ic  =  3IQ 


and  add  and  subtract 
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3la  3Ib 

L13  3t~  +  L13  W  from  eq*    (A"la) 

aia  aib  3i 

Ll3  at"  +  L13  W  +  L13  at-     from  eq-  (A"lb) 
aib  31 

and  L13  at-  +  L13  3t~~  ^rom  eq*   (A~ic)* 

The  result  is 

3V*  3Ia  3Ih 

"  JT  -  RUla  +   <LirL13>   3T  +   <L12-L13>   3T  +  ^O  + 

3(L00+L13>  W-  (A"3a) 

3Vb  3Ia  3Ib  31 

"  3z~  =   (L12"L13^   3t~  +  R22Jb  +   (L22"L13^  W  +   ^L23'L13)   3t~  + 

3FWo  +  3(Loo+Li3>  Ir        <A-3b> 


3V  31  31 


c 


=   (Lo,-L,J  5^+  R«I,  +   (L„+L,o)  t^+  3RnnI 


3z         vu23  u13'   at         "331c       vu33  L13y   at         jrW00 

3(L00+L13>  ar  (A"3c> 

Equations   (A-2)  and   (A-3)  are  the  relations  describing  the  "balanced 
mutual"  model   of  Figure  2.2  with: 

Lll        Lll        L13  L12       L12       L13  L10       L10 

L'=L-L  L'=L-I  C*      =   C 

L22       L22       L13  L23       L23       L13  ^20       L20 

l'=l  -I  I     '  3     I  +1  C  =      C 

L33   L33   L13         L00   L00   L13         L30   L30 

C12  =  C12  Rll  =  Rll  R00  =  R00 

C13  =  C13  R22  =  R22 

C23  =  C23  R33  =  R33  ^A"4^ 

A. 3     Transformation  of  the  General   Non-Transposed  Model   to  a  Transposed 
Model 

One  transposition  cycle  of  a  transposed  three-phase  transmission 

line  can  be  visualized  physically  as  shown  in  Figure  A.l.     In  the 
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transposed  line  let: 

L.j  =  self  inductance  of  a  phase  =  Ls 

L-.  =  mutual  inductance  between  two  phases  =  L.. 

C.  =  capacitance  to  ground  of  a  phase  =  CQ 

C .  .  =  mutual  capacitance  between  two  phases  =  cL 

R..  =  self  resistance  of  a  phase  =  Rs 

Then  from  figure  (A-l)  it  is  seen  that  in  terms  of  the  parameters  of 

Figure  2.1 : 

,T         ,T         .T      ,  Lll   *  L22  *  L33       ,T  ,.   _   » 

Lll    '     22  ""     33  "  3  "'  LS  lA-aa; 

iT      -iT         iT      .  L12  *  L13  +  L23  ,.  ,T  ,.   ,.. 

12  "     13         23  3  "     M  IA-3DJ 

rT         rT         rT         C10  *  C2Q  *  C30       rT  ,.   ,   v 

L10  "     20  "'  L30  "  3  "  L0  lA-sc; 

c     +  c     +  c 
rJ     -  cJ      ■  rT     -     12       L13  +  L23  ..  fT  ,.       * 

12  "     13  "     23  "  3  "'     M  IM-oa; 

R!l   ■  R?2  ■  R33  "  Rl1   +*2l     +  "33  ■  RS  <A"5e> 

The  significance  of  these  relations  can  be  seen  more  clearly  by  modeling 

a  transposition  cycle  of  the  transposed  line  with  the  model  of  Figure 

2.1.     The  differential   equations  corresponding  to  equations   (A-l)  with 

I     +  I.    +  I     =  3In  are: 
a        b        c  l) 

9Va         T  31  31  31  31 

"  IT  =  Rsh  +  LsV+  LS  runaT+  ^oo'o  +  3Loo  aT  (A"6a> 

3V.  T  3I  =         T  T  3Ih         T  31.  3In 

"  JT  -  4  3T  +  Rs'b  +  Ls  3T  +  4  3T  +  3Wo  +  3Loo  ST         (A"6b) 

3V  T  31  31  T  T  31  31- 

"  3T  ■  4  3T  +  4  5T  +  Rs'c  +  4  5T  +  ^oo'o  +  34o  3T  •       <A-6c> 
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T  3la  b 
Adding  and  subtracting  L,  '. '  as  appropriate  in  equations  A-6  obtain: 

3Va    T      T    T  9Ia  T       3IQ 

-  JT  -   RSJa  +  <4  -   W  JT  +  *RQQl0   +  3(LM  +  L00>  W-  (^ 

3V,  T  T         T     31.  T  3In 

"  JT  -  RS!b  +    <LS   "  LH}    3T  +  3fWo  +  3K  +  L00>  JT  <A"7b> 

3Vr         T  T         T     3I_  T  3In 

"  3T  "  Rs!c  +  <Ls  "  hi)  JT  +  3Wo  +  3(lm  +  Loo>  JT  <A"7c> 

In  the  balanced  steady-state  IQ  =  0;  the  currents  in  the  transposed  line 
are  positive  sequence  only,  and  equations   (A-7)  reduce  still   further  to 

T       T 

three  decoupled  "single-phase  equivalent"  equations.  With  Ls-LM=Lp  s  = 
positive-sequence  inductance  and  Rs  =  Rp  s  =  positive-sequence  resistance 

3V  31 

3V,  31. 

-ar"Rp.s.Ib  +  Lp.s.  w-  <A"8b> 

3Vr  31 

-JT=RP.s.lb  +  lP.s.W-'  (^ 

If  only  zero-sequence  currents  flow  in  the  line,   I     =  I,    =  I    ,  and 
equations    (A-7)  reduce  to 

3Va           T  T  T                   31 

"  JT  -  (RS  +  Va  +   (4  +  2L!1  +  3L00>  JT                     <A"9a> 

3V,             T  T  T                    31. 

"   3^  "  <RS  +  Vb  +   (LS  +  2LM  +  3L00>  JT                      (A"9b) 

"  JT  =  (Rs  +  3Roo)rc  +  (LS  +  K  +  3L00>  JT  ^9c) 

Equations    (A-8)  and   (A-9)  describe   (in  terms  of  The  Voltage  Law)  the 
symmetrical   components  equivalent  circuits  of  the  transposed  line 
described  by  equations    (A-7).     The  transposed  line  corresponding  to 
equations   (A-7)   is  shown  in  Figure  2.3. 
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The  Current  Law  equations  for  the  transposed  line  corresponding 
to  equations  (A-2)  can  be  seen  by  inspection  of  Figures  2.1  and  A-l : 

dK       t     t  9V,    T  9Vh    t  dyr 
31.      T  3Va     T     T  3V,    T  3V 

"  5T  ■  -  ch  3T  +  <  +  2Cn>  W-  -   cn  3T  <A"10b) 

31       T  3V     T  3V,     T     T  3V 

-  w-  -  -  ca  ir  -  cn  w- +  (co  +  2CJ>  w-  <A-10c» 

T  9Va  b 
Setting  V  +  V.  +  V  =  3Vq  and  adding  and  subtracting  C,,  ~^-  from 

equations  (A-10)  as  appropriate: 

-  ir  ■  ^ +  3CS>  ar  -  3CS  ?r  (A-lla) 

31.     T     T  3V.     T  3V 

31      T     T  3Vr     T   3Vn 

-JT-  (c5  +  39  3T-  3CM3-r  <A-11c> 

In  the  balanced  steady-state  V0  =  0;   the  voltages  are  positive  sequence 

only,  and  equations   (A-l 1 )   reduce  to  three  decoupled  "single-phase 

equivalent"  equations.     With   (Cq  +  3C.,)  =  Cp  s     =  positive  sequence 

capacitance,  equations   (A-l 1 )  become: 

31  3V 

=  Cn  c     tt^-  (A-12a) 


a 


3z  P.S.    3t 

31,  3V. 

31  3Vr 

-JT=CP.S.  3T  <A-12c) 

If  only  zero-sequence  voltages  are  applied  to  the  transposed  line, 

V  =  V.  =  V  .  And  with  Cn  =  C-,  c  =  zero  sequence  capacitance, 
a    b    c  0    Z.S.        ^       r 

equations  (A- 11)  become: 
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9Ia 

31 


3z 

31 


9Va 

z. 

,s. 

3t 
9Vb 

z, 

,s. 

3t 
3Vc 

1  =  C7  c  ^  (A-13b) 


1  =  C7  c  ^  (A-13c) 


3z     Z.S.  3t 

Equations  (A-12)  and  (A-13)  describe  (in  terms  of  The  Current  Law)  the 
symmetrical  components  equivalent  circuits  of  the  transposed  line 
described  by  equations  (A-ll).  The  transposed  line  model  corresponding 
to  equations  (A-ll)  is  shown  in  Figure  2.3. 
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APPENDIX  B 

MODEL  IMPEDANCE  DERIVATIONS 

B.1.1     Derivation  of  L-Section  Model   Impedance- Inductor  Input- 
Shorted  Termination  ~~~  

Using  the  m  section  line  shown  in  Figure  3.3,  write  difference 
equations  near  the  nth  section   (0<n<m): 

Vn  "  Vn-1  "  Ls  3T  +  Rs'n  (B-D 

dV 
Vl   "  'n  =  Cs  dT  (B-2) 

Vl   -'n^sJT1*  Vn+1  (B-3) 

Subtracting   (B-l)  from  (B-3): 

Vl  "  2Vn  +  Vl  - L  at  <VrV +  MWV  •        (B-4) 

Substitute  (B-2)  into  (B-4): 

d2V„  dV 

Vn-1  "  2Vn  +  Vl  "  LsCs  TT  +  RsCs  iT  ■  (B-5) 

dt 

Assuming  solutions  of  the  form  V     =  Re[V     est]: 

n  n 

Vl   -  Z\  +  Vl  ■  *\V„  +  s  RsCs5n  <B-6) 

*  ±nA 

Next  assume  Vn  =  K+    e"       (where  X  may  be  a  complex  constant, 

*=*r+j*.j)  and  substitute  into  equation  (B-6): 

e±(n-1)X  .  2  e±nX  +  e±(n+l)X  .(  2  +  )  ±„X 

s  s  s  s' 

e~A  -  2  +  eA  =  s2L  C     +  s  R  C 

s  s  s  s 

s2|-sCs       s  RsCs 
cosh  X  =  1  +  1-2-  +        s  s  (b-7) 
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Therefore,  provided  A  satisfies   (B-7),  the  solutions  to  equation  (B-5) 
take  either  of  the  forms: 

Vn  =  Re[Vn  est]  =  Re[(K+  enA+  K_  e"nA)est]         (B-8) 

Vn  =  Re[(A  cosh  nA  +  B  sinh  nA)eSt]  (B-9) 

To  obtain  the  current,   I   ,  use  equation  (B-l)  assuming  solutions  of  the 


form: 


In  =  Re[In  est]  (B-10) 

Vn  =  Re[Vn  est]  .  (B-ll) 


Then  (B-l )  becomes: 

Vs  Ls+Rs>  ■  Vl   "  \ 
V     -  V     , 

'H - *"+  s ";  •  <B-,2» 

s  s 

In  terms  of  equation   (B-9)  equation   (B-12)  can  be  written  as: 

T         A(cosh  nA  -  cosh(n-l)A)  +  B(sinh  nA  -  sinh(n-l)A)  /D  loX 

ln=  Rs  +  s  Ls (B"13) 

If  the  short  is  connected  at  the  zeroth  node,  then  VQ  =  0  and  from 

equation   (B-9)  it  is  seen  that  A  =  0,  so  that  elsewhere  on  the  line: 

V     =  B  sinh  nA   .  (B-14) 

n 

And  from  equation   (B-l 3): 

C         B(sinh  nA  -  sinh(n-1 )A)  /D  1C\ 

I     =  r     +  s  L  *  lB"lbj 

s  s 

Finally,  the  impedance  looking  into  the  mth  node  -  or  driving  point 

terminals  -  is  seen  to  be: 

(R     +  s  L   )sinh  mA 
Zm(s)  =  sinh  mA  -  sinh(m-l)A  (B"16) 

Where  it  is  noted  again  that  A  must  satisfy 
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s2LC         s   R.C 


cosh  A  =  1   +  — sM-  +  !-§. 


(B-7) 


In  the  case  that  the  model   is  lossless   (Rs  =  0),  equations   (B-16)  and 
(B-7)  reduce  to: 


jcoL     sinh  mA 


Zm(ja,)  =  sinh  mA  -  sinh(m-l)A  (B-17) 

u2lsCs 
cosh  A  =  1 p.  .  (B.18) 

The  solution  for  the  poles  of  (B-17)  is  indicated  in  Appendices  B.1.2, 
C.2.1,  and  D.l. 

B.1.2     Solution  for  the  Poles  of  the  L-Section  Model-Inductor  Input- 
Shorted  Termination  ~~  " ~~~ —  — 

The  solution  for  the  poles  of  equation  (B-16)  or  (B-17)  rests 
upon  recognition  of  the  fact  that  the  denominator  of  (B-16)  and  (B-17) 
is  a  polynomial   in  cosh  A  of  degree  (m-1)  for  all  m.     That  this  is 
true  can  be  seen  by  expansion  of  the  denominator  for  the  first  few 
values  of  m.     The  expansions  of  sinh  mA  up  to  m  =  7  are  shown  in 
Table  B.l. 

Table  B.l 
Expansions  of  sinh  mA 

sinh  A     =  sinh  A 

sinh  2A  =  sinh  A(2  cosh  A) 

sinh  3A  =  sinh  A(4  cosh2A-l) 

sinh  4A  =  sinh  A(8  cosh3A-4  cosh  A) 

sinh  5A  =  sinh  A(16  cosh4A-12  cosh2A+l ) 

sinh  6A  =  sinh  A(32  cosh5A-32  cosh3A+6  cosh  A) 

sinh  7A  =  sinhA(64  cosh6A-80  cosh4A+24  cosh2A-l) 
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The  expansions  listed  in  the  table  were  obtained  by  repeated 
application  of  the  identity  sinh  mX=sinh(m-l )X  cosh  X+cosh(m-l)X  sinh  X. 
The  use  of  this  identity  to  obtain  sinh  mX  requires  the  preceeding 
recursion  terms  sinh(m-l)X  and  cosh(m-l)X.  The  term  sinh(m-l)X  is 
simply  the  preceeding  term  in  the  table,  but  cosh(m-l)X  must  be  obtained 
from  a  separate  process  of  expansions  of  cosh  mX.  Table  B.2  shows  the 
first  few  expansions  of  cosh  mX  generated  in  the  loop  required  to  get 
sinh  mX. 

Table  B.2 
Expansions  of  cosh  mX 
cosh  X  =  cosh  X 

cosh  2X  =  2  cosh2X-l 

3 
cosh  3X  =  4  cosh  X-3  cosh  X 

cosh  4X  =  8  cosh4X-8  cosh2X+l 

cosh  5X  =  16  cosh5X-20  cosh3X+5  cosh  X 

cosh  6X  =  32  cosh6X-48  cosh4X+18  cosh2X-l 

cosh  7X  =  64  cosh7X-112  cosh5X+56  cosh3X-7  cosh  X 


Using  Table  B.l,  the  first  seven  polynomials  set  equal  to  zero  in  the 

denominator  of  equations  (B-16)  and  (B-17)  are  listed  in  Table  B.3. 

Table  B.3 
Expansions  of  Impedance  Function  Denominator 

0  =  1-  cosh  X  =  0  (m=l ) 

(m-2) 
(m=3) 

0  =  '  -  4  <  sh  X-4  <  osh^X+S  cosh~X     etc. 


0  = 
0  = 


-  cosh  X  =  0 

-  2  cosh  X 
2 

+  2  cosh  X-4  cosh  X 

9        3 

-  4  cosh  X-4  cosh^X+S  cosh  X 

+  4  cosh  X-12  cosh2-8  cosh3X+16  cosh4X 

-  6  cosh  X-12  cosh2X+32  cosh3X+16  cosh4X-32  cosh5X 
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The  equations  listed  in  Table  B.3  are  the  equations  which  must  be  solved 
to  obtain  the  poles  of  L-section  models  of  from  one  to  seven  sections. 
In  obtaining  the  equations  of  Table  B.3  the  sinh  X  term  is  divided  out 
since  sinh  X  t   0.  A  more  thorough  discussion  of  the  constant  X  is  now 
in  order.  In  general,  X  ■  X  +  jX . ,  so  that  in  equation  (B-18),  for 


example, 


w  LsCs 

cosh  X  ■  cosh  X„,  cos  X.  +  i  sinh  A  sin  X.  =  1 * .   (B-19) 

r  i  r  l  c 


Equation   (B-19)   implies: 


<A_ 


cosh  \     cos  X.   =  1 %—-  (B-20) 

r  l  c 


sinh  Xp  sin  Xi   =  0  (B-21) 


But,  again, 


2, 

03  L   C 


s  s 


cosh  X  =  cosh  Xr  cos  X.   =  1 ^-  (B-23) 

and  it  is  seen  that  cosh  X  is  real,  and  X  is  pure  imaginary  (=  jX.)  if 

oj2L  C 

\3-  >   0.  With  X  =  jXif  Xp  =  0  so  that  finally, 

oj2L  C 
cosh  X  =  cosh  jXi  =  cos  Xi  =  1 2^"  •  (B-24) 

Equation  (B-24)  shows  that  all  values  of  cosh  X  must  be  real  and  lie 
between  minus  one  and  plus  one.  This,  in  turn  implies  that  all 
solutions  to  the  "equation"  in  the  denominator  of  (B-17)  must  also  be 
real  and  lie  between  minus  one  and  plus  one  i.e., 

-1.0  <  [cosh  X]root  <  1.0  .  (B-25) 

But  the  denominator  of  equation  (B-17)  is  the  same  as  that  of 
equation  (B-16),  so  equation  (B-25)  applies  for  the  lossy  line  also. 
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Thus,  the  equations  listed  in  Table  B.3  are  polynomial   equations  in  the 
real   variable  cosh  A  for  which  there  are  numerous   "root-finding"  sub- 
routines in  any  good  computation  facility.     With  the  roots  of  equations 
like  those  in  Table  B.3  readily  available,  equations   (B-7)  or  (B-16) 
can  be  solved  for  s  or  u  respectively  to  give  the  poles: 


ss.c.   ■  -  ?L  ±V(fr)2  "  J^  H-tcosh  A)root]  (B-26) 

hi  '   1,2,3 


0)- 

s 


.c.  ■  i  Jb  ^-(cosh  *  W  <B-27> 

m  =   1,2,3 

Where  use  has  been  made  of  the  fact  that: 

Rc  =  %r  »  Lc  =  ^r  »  and  cc  =  ^r  (B-28) 

s   m    s   m       s   m 

Equations  like  those  in  Table  B.3  are  seen  to  be  solved  easily  enough, 

but  the  problem  of  how  to  generate  the  equations  up  to  fifty  or  more 

sections  (m)  remains.  The  details  of  the  generation  problem  are  left 

to  Appendix  C. 

B.2     Derivation  of  L-Section  Model    Impedance,  Capacitor  Input,  Shorted 
Termination 

Using  the  m-section  line  shown  in  Figure  3.6,  write  difference 

equations  near  the  nth  section   (0  <  n  <  m): 

Vl  -  Vn  =  Ls  dT1  (B"29) 

dIn     1 

n         n-1         s     dt 

dV 

I     -   I     t   =  C     -,=£■  (B-31) 

n         n-1         s  dt 
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Subtracting   (B-30)  from  (B-29): 

Vn+1  "2Vn  +  Vl  =  Ls3t  "n"  W'  (B-32) 

Substitute   (B-31)   into   (B-32): 

d2V 


.fl  -  2Vn  +  Vl  "  LsCs  -r 


Similar  to  tiie  development  in  Appendix  B.l.l  assume, 

Vn  =  Ke[Vn  e±ja)t] 


Vn  =  K±  e±n* 

nd  obtain: 

2,    r 

d)L  C 

cosh  A  =  1   -         s  s 

(B-33) 

Vn  =  A  cosh  nA  +  B  sinh  A  (B-34) 

To  obtain  the  current,  In>  assume  Ip  =  Re[In  ejwt]  in  equation  (B-29): 

jwLs*n  =  Vi  -  Vn  (B.35) 

1 
Jn  =  JJT  tA  cosh(n+l)A  +  B  sinh(n+l)A  -  A  cosh  nA  -  B  sinh  nA] 

( R-  "3fi  ^ 
With  a  short  at  the  zeroth  node  VQ  =  0  implies  A  =  0,  and 

Vn  =  B  sinh  nA  (B_37) 

B 
Jn  =  M~  Lsi nh(n-i-l  )A  -  sinh  nA].  (B-38) 

s 

Finally,  for  the  driving  point  impedance  at  the  mth  node, 

7  ..,  .      sinh  mA  ,    _  „  , 

Zs.c.   -  JU)Ls  sinh(m+l)A-sinh  mA       m  =  1*2'3 <B"39) 

and  using  equation  (B-33): 
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us.c.  =iyfetl-(cosh  A)roQt]   m-  1,2,3. 


(B-40) 


The  solution  for  the  poles  of  equation  (B-39)  is  indicated  in  Appendices 

C.2.2  and  D.2. 

B.3     Derivation  of  L-Section  Model    Impedance,  Open  Termination 

Using  the  m  section  line  shown  in  Figure  3.3,  but  with  the  short 
removed,  R     =  0,  and  s  =  ju>,   the  derivation  proceeds  as  in  Appendix 


B.l .1   to  obtain: 


V     =  A  cosh  nX  +  B  sinh  nX 


1 


n  "  jooL 


(B-41) 


A[cosh(n-l)X-cosh  nX]  +  B[sinh(n-1  )X-sinh  nX]  (B-42) 


But  now,   Iq  =  0,  and  equation   (B-42)   implies, 


P  _  cosh  X-l 


sinh  X 
Substituting   (B-43)   into   (B-41)   and   (B-42): 


(B-43) 


Vn  =  A 


"sinh(n+l  )X-sinh  nX 


A 


n  "  jwL 


sinhX 
2  sinh  nX-sinh(n+l )X  -  sinh(n-l)X" 


sinh  X 


(B-44) 
(B-45) 


So  that  looking  into  the  mth  section: 


Zo.c.  ■  J"Ls 


As  in  Appendix  B.l .1 


["        sinh(m+1  )X-sinh  mX 
[2  sinh  mX-sinh(m+l)X-sii 


nh(m-l)X_ 


cosh  X  =  1 


2i    r 
o)LsCs 


from  which 


Vc.  -lACl   "   (GOShX)root 


(B-46) 


(B-18) 


(B-47) 
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The  solution  for  the  poles  of  equation   (B-46)  is  indicated  in  Appendices 
C.2.3  and  D.3. 

B.4     Derivation  of  TT-Section  Model   Impedance,  Shorted  Termination 
Using  Figure  B.l,  write  difference  equations  near  the  nth 
section: 


dl, 

n       Vl       Ls  dt 

dvn 

I     _  T     =  r    D. 

n       L*      Ls  dT 


I*  -   I   .,   =  C 


dV 


n+1 


*       An+1       °s  ~dt 
Substitute  (B-50)  into  (B-48)  to  eliminate  I*: 

dV 


(B-48) 
(B-49) 

(B-50) 


V„  -  V  „  =  L 


C    n+1  i  T 
Ls  dt    'n+1 


n    n+1    s  dt 
Substitute  (B-50)  into  (B-49)  to  eliminate  I*: 


(B-51) 


dV 


I  -  C. 


n+1 


dV. 


I  .-,  =  C 


(B-52) 


n    s  dt    An+1   "s  dt  * 
Write  the  equivalent  to  equation  (B-51)  for  the  next  section  toward  the 
left: 


V  ,  -  v  =i  — 
n-1   vn   Ls  dt 


dV 

'  — +  I 
's  dt    ln 


(B-53) 


From  equation  (B-52)  write; 


'n+1  "  <n  "  Cs  ft  [Vn+1  +  V  • 
Substitute  (B-54)  into  (B-51)  to  eliminate  I 

dl 


d2V 


n+1 
2 


V.  -  V_,  =  L.C 


n    n+1    s  s   dt 
Rewrite  equation  (B-53): 

V_  =  L.C. 


+  L_rn.L_C.  <^(v    +VJ. 


•s  dt    usus  dt2  v  n+ 


d2V 


dl. 


+  L 


n-1    n   us  s  d  2    Ls  dt 


(B-54) 


(B-55) 


(B-56) 
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dl 
Substitute  (B-56)  into  (B-55)  to  eliminate  l-r^-   • 

dt 

d2V 


Vl  "  2Vn  +  Vn-1  =  2LsS  -^  (B-57) 

dt 


Solve  equation  (B-57)  as  in  earlier  cases  to  obtain: 


cosh  A  =  1  -  W2L  C 


1 


and 


(B-58) 


and 

Vn  =  A  cosh  nA  +  B  sinh  nA  .  (B-59) 

To  find  the  current  use  equation  (B-56)  assuming  solutions 

I"  =  Re[in  ejwt] 

Vn  =  Re[Vn  ejwt] 

to  obtain: 

'n  "  JZ  tf„-l  "  \   +  "VsV  ■  (B-60) 

Expanding  equation  (B-60): 


=  '  j^L  ^A  sinh  X  sinh  nA  +  B  cosh  nA  sinn  *L      (B-61) 


With  a  short  at  the  zeroth  node,  VQ  =  0  implies  A  =  0  so  that 

Vn  =  B  sinh  nA  (B-62) 


7     -       -i\>i  sinh  nA  . 

n  ~  ■  JwLs  cosh  nA  sinh  A   '  (B"63) 

At  the  mth  node  or  driving  point  terminals  : 

7     _        •  ,  sinh  mA  .     „.. 

m  "  "  JwLs  cosh  mA  sinh  A     '  <B"64) 

The  solution  of  equation   (B-64)  for  the  poles  is  discussed  in 

Appendices  C.2.4  and  D.4.     The  results  of  the  solution  are  the  values 

of  [cosh  A]roQt  which  give  the  poles  of  (B-64).     Also,  using 
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equation   (B-58), 


or 


1   -   (cosh  X)       . 

»s.c.  -v/ LC  (B"65> 

s  s 

m    /2L1  -  (cosh  X)rootJ 

us.c.  =  W cc '  (B"66> 

Equation  (B-66)  follows  from  (B-65)  since 

Ls-r    and    cs*?sr  •  <B-67> 

B.5     Derivation  of  ir-Section  Model    Impedance  Open  Termination 

The  derivation  of  the  open  circuit  case  is  the  same  as  that  of 
the  short  circuit  case  through  equation  (B-61).     Now,  at  the  "0"th  node, 
IQ  =  0,  and  equation   (B-61)  implies  B  =  0.     Thus 

I     =  -  tV  [A  sinh  X  sinh  nX]  (B-68) 

Vp  =  A  cosh  nX  (B-69) 

V 
7  _  m  _   .  -     cosh  mX  ,D  -,n\ 

Zm  -  j-  -  -  juL  sinh  x  s1nh  mX  (B-70) 

m 


Again, 


m    /2L1  -  (cosh  XT~ 
m    /    L  'root 


Jo.c.       £v  LC 

The  solution  of  equations   (B-70),   (B-71),  and  (3-1    )  is  discussed  in 
Appendices  C.2.5  and  D.5. 


(B-71) 
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APPENDIX  C 
SOLUTIONS  FOR  THE  POLES  OF  THE  MODEL  IMPEDANCE  FUNCTIONS 
C.I  Model  Impedance  Functions 

The  model  impedance  functions  obtained  in  Appendix  B  for  which 
solutions  are  sought  are  shown  in  Table  C.I.  For  the  purposes  of 
Appendix  C,  "solutions"  refers  to  those  values  of  cosh  A  in  the 
denominator  of  the  impedance  function  which  result  in  the  poles  of  the 
function.   In  Appendix  B.1.2  it  was  shown  that  the  denominator  of 
equation  (B-17)  is  a  polynomial  in  the  real  variable  cosh  A.  In 
addition  it  was  shown  that  the  polynomial  had  real  coefficients  and 
that  all  solutions  for  [cosh  *]root  should  be  such  that 

-1.0  <  [cosh  A]roQt  <  1.0  . 

Since  the  denominator  of  equation  (B-17)  is  a  polynomial  in  cosh  A, 
it  follows  that  the  denominator  of  all  the  impedance  functions  is  a 
polynomial  in  cosh  A  with  all  the  properties  of  the  polynomial  in 
equation  (B-17).   It  is  the  denominator  polynomial  set  equal  to  zero 
which  results  in  the  poles  of  the  respective  impedance  function.  If 
the  denominator  polynomials  [D(cosh  A)]  can  be  generated  for  any 
number  of  sections,  m,  then  it  is  a  simple  matter  to  call  upon  root 
routines  in  the  computer  library  to  solve  for  the  roots  of  the  denominator 
equation,  D  =  0.  An  inspection  of  Table  C.l  and  reference  to  Appendix 
B.1.2  shows  that  the  "generation"  of  D  involves  first  the  generation  of 
sinh  mA  and  cosh  mA  for  any  m.  In  Appendix  B.1.2  it  was  shown  that 

sinh  mA  =  sinh  A[Polynomial  in  cosh  A]  (C-l) 


141 


Model 


Driving  Pt.  Impedance  Function 


m 

= 

jcoL     sinh  mX 

L-Section 
Inductor  Input 

sinh  mX-sinh(m-' 

nx 

Shorted  Termination 

ijb  [1-(cosh 

X)root] 

CO 

= 

jcoL     sinh  mX 
s 

L-Section 
Capacitor  Input 
Shorted  Termination 

sinh(m+l )  -sinh 

m 

ijk  ^-(cosh 

A)root^ 

Zm 

m 

CO 

= 

jcoL  [sinh(m+l)X-sinh  X] 

L-Section 

Open  Termination 

2  sinh  mX-sinh(m+l )X-sinh(m-l )X 

ijhv-^ 

XVoot^ 

Zm 

m 

CO 

= 

-jcoL     sinh  mX 

7r-Section 

Shorted  Termination 

cosh  mX  sinh  X 

KL^[Mcosh 

AVoot] 

Zm 

m 

CO 

= 

-jcoL     cosh  mX 

Tr-Section 

Open  Termination 

sinh  X  sinh  mX 

l7lC^-(cosh 

XVoot^ 

:b-i7) 

;b-27) 

;b-39) 
;b-4o) 

;b-46) 
;b-47) 

;b-64) 
;b-66) 

;b-7o) 

:b-7d 


Table  of  Impedance  Functions  for  Which 
Poles  are  Sought 

Table  C.l 
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and 

cosh  mX  =  [Polynomial  in  cosh  X]  .  (C-2) 

Denote  the  polynomial  in  equation  (C-l)  by  B  =  B(cosh  X)  and  that  in 
equation  (C-2)  by  A  =  A(cosh  X)  so  that 

sinh  mX  =  sinh  X  Bm(cosh  X)  (C-3) 

cosh  mX  =  Am(cosh  X).  (C-4) 

Then  it  follows  that  in  all  the  impedance  functions  of  Table  C.l  the 
denominator  has  the  form 

D(cosh  X)  =  Linear  Combination  of  A(cosh  X)  and  B(cosh  X). 

(C-5) 

A  routine  must  be  found  which  will  enable  the  computer  to  generate  A 

and  B  .  combine  them  to  form  D  =0,  and  then  solve  for  the  values  of 

m  m 

[cosh  Xl   .  for  use  in 
root 


w  =  I  Jb  [1  "  (cosh  X)root]  *  (C_6) 

C.2  Generation  of  Impedance  Function  Denominator  Polynomial 
In  general 

sinh  mX  =  sinh(m-l)X  cosh  X  +  cosh(m-l)X  sinh  X, 

so  that  in  the  terminology  of  Appendix  C.l, 

sinh  X  Bm  =  sinh  X  Bm  ,  cosh  X  +  Am  ,  sinh  X, 
m         m- 1         m- 1 


and 


Bm  =  Bm  ,  cosh  X  +  A  ,  .  (C-7) 

m    m-1         m- 1 


Similarly 

cosh  mX  =  cosh(m-l)X  cosh  X  +  sinh(m-l)X  sinh  X, 
so  that 
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Am  =  K,   i  cosh  A  +  Bm  ,  sinh2A, 
m   m- 1         m- 1 


or 


Am  =  Am-1  cosh  X  +  (cosh2A-l)Bm_1  .  (C-8) 

Throughout  the  remainder  of  this  appendix  denote 

cosh  A  by  "C" 
and 

sinh  A  by  "S" 
so  that  equations  (C-7)  and  (C-8)  become: 

Bm  =  C  Bm  1   +  Am  1  111  =  1,2....         (C-9) 

m     m-1    m-1 

Am  =  c  Am  i  +  (c2~1)Bm  ,    m  =  1,2....     (C-10) 
m     m-1       '  m-1 

With  m  =  1,2 ,  comparison  of  equations  (C-9)  and  (C-10)  with  Tables 

B.l  and  B.2  indicates  that  with  AQ  =  1  and  BQ  =  0  equations  (C-9)  and 

(C-10)  form  a  recursion  loop  for  generating  all  B  and  A  .  In 

equations  (C-9)  and  (C-10)  it  is  the  coefficients  of  the  powers  of  cosh  A 

which  are  needed  to  solve  for  the  roots  of  the  "D"  polynomial,  so  that 

in  B  and  A  it  is  the  coefficients  which  are  of  interest. 

m     m        

In  their  most  general  form  the  polynomials  A(m)  and  B(m)  take 
the  form  (for  m  =  1,2 ) : 

A(m)  =  am+1(m)Cm  +  am(m)Cm_1  +  *} (m)  (C-ll) 

B(m)  =  bm(m)Cm_1  +  bm_] (m)Cm"2+ b] (m)  .         (C-12) 


Where 


m  =  number  of  sections  in  the  model 

C  =  cosh  A 

A(m)  =  "A"  polynomial  when  model  has  m  sections 

B(m)  =  "B"  polynomial  when  model  has  m  sections 
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a. (m)  =  coefficient  of  C   when  model  has  m  sections 
b..(m)  =  coefficient  of  C1"  when  model  has  m  sections 
The  subscripts  of  the  coefficients  begin  at  "1"  in  order  to  be  compatible 
with  the  subscripting  convention  in  the  Fortran  IV  language.  Note  from 
Tables  B.l  and  B.2  that 

B(l)  =  1  implies  b^l)  =  1, 

A(l)  =  C  implies  a2(l)  =  1,  (C-13) 

and    a, (1)  =  0. 

To  compute  B2,  for  example,  use  equations  (C-9),  (C-10),  and  (C-ll): 

B2  =  B(2)  =  C  B(l)  +  A(l) 

B(2)  =  CCb^l)]  +  [a2(l)  C  +  a^l)] 
B(2)  =   [a2(l)  +  b^l)]  C  +  a^l) 


Where 


and 


In  general : 


[a2(l)  +  b^l)]  =  b2(2) 


a^D  =  ^(2) 


B(m)  =  C  B(m-l)  +  A(m-1 ) 

=  C[bm.1  (m-1  )Cm_2+bm_2(m-l  )Cm~3. . .  .b1  (id-1  )] 

+  [am(m-l)Cm'1+an  , (m-1 )Cm"2 a, (m-1)] 

m  iii- 1  i 


Combining  coefficients  of  like  powers: 

.m-1 
,m-2 


B(m)  =  +  [bm_.,(m-l)  +  am(m-l)]Cn 


+  [bm.2(m-l)  +  Vl(m"1)tf 
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+  [b^m-1)  +  a2(m-l)]C 

+  [a^m-1)]  .  (C-14) 

Inspection  of  equation  (C-14)  shows  that 

b^m)  =  [b^Cm-l)  +  a^m-1)]   2  <_  i  <  m        (C-15) 
b^m)  =  a^m-1)  .  i  =  1  (C-16) 

A  similar  procedure  applied  for  A(m)  gives: 

a^m)  =  [a^^m-1)  +  bi_2(m-l)]         i=m,  m+1  (C-17) 

a^m)  =  [a-^dn-l)  +  b1-2(m-l  )-b.  (m-1 )]  i=2,3,. . .  ,m-l 

(C-18) 

a2(m)  =  [a^m-1)  -  b2(m-l)]  i=2     (C-19) 

a^m)  =  L-b^m-1)]  .  i  =  l      (C-20) 

Given  the  initial  values  for  A(l)  and  B(l),  equations  (C-15)  through 
(C-20)  can  be  implemented  in  Fortran  IV,  and  the  computer  made  to 
generate  and  store  the  coefficients  of  the  "A"  and  "B"  polynomials. 
Having  the  coefficients  of  "A"  and  "B"  up  to,  say,  fifty  sections 
(m=50),  the  "0"  polynomial  can  be  generated  by  simple  algebra. 
C.2.1  Polynomial  for  L-Section  Model,  Inductor  Input,  Shorted 

For  the  case  of  the  L-Section,  inductor  input  model  with  a 
shorted  termination,  the  equation  for  which  the  roots  are  sought  is 
(from  equation  B-17) 

sinh  mX  -  sinh(m-l )X  =  0 
sinh  X[B(m)  -  B(m-l)]  =  0. 
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The  case  A  =  0  has  no  meaning  for  this  model  and  so 

B(m)  -  B(m-l)  =  0,  (C-21  ) 

and  equation  (C-21)  is  seen  to  be  the  relation  which  gives  the  poles  of 
equation  (B-17).  Thus 

U(m)  =  B(m)  -  B(m-l)  =  0.  (C-22) 

Equation  (C-22)  must  also  be  written  for  implementation  into  Fortran  IV 
in  the  same  way  that  A(m)  and  B(m)  were  implemented.  Thus: 

D(m)  =  dm(m)Cm_1  +  d|n_1(m)Cm"2 d£(m)C  +  d^m) 

=  B(m)  -  B(m-l) 

=  ^(mjc"1"1  +  bm_1(m)Cm-2  +....b2(m)C  +  b^m)] 

-  [bm_.,(m-l)Cm-2  +  Bm_2(m-l)Cm-3  +....b2(m-l)  +  b^m-1)]. 

Combining  coefficients  of  like  powers: 

Jn-1       ru       /_>       L       /     ,  M^m-2 


U(m)  =  bm(m)Cm-'   +  [b^m)  -  bm_1  (m-1  )]Cm^  + 
[bm.2(m)  -  bm_2(m-l)]Cm-3  + 


,..[b2(m)  -  b2(m-l)]C  +  [b]  (m)  -  b^m-1)]. 
So  that  identifying  coefficients  of  D(m): 

d.(m)  =  b.(m)  -  b.(m-l)  i   =  1 m-1  (C-23) 

d   (m)  =  b   (m)    .  (C-24) 

Thus  the  computer  can  be  made  to  generate  the  coefficients  of  the  "D" 
polynomial.  The  form  of  the  coefficients  given  in  equations  (C-24)  and 
(C-25)  are  explicitly  suited  for  several  "root"  subroutines  available 
in  the  IBi-1-360  Scientific  Subroutine  Package  (SSP)  at  the  Computation 
Center  of  M.I.T.  Equations  (C-15)  through  (C-20)  and  (C-23)  and  (C-24) 
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were  implemented  in  Fortran  and  solved  with  the  IBM- 360  as  shown  in 
Appendix  L).l.      In  Appendix  C.l    the  roots  of  D  are  put  into  equation 
(B-26)   to  obtain  the  poles  of  a  100  mile  long  typical    765  kV  line. 
The  poles  of  the  model    thus  obtained  are  compared  to  the  poles  of  the 
actual    line  to  obtain  the  pole  error  plots  of  Chapter  3. 
C.2.2     Polynomial   for  L-Section  Model,  Capacitor  Input,  Shorted 

For  the  case  of  the  L-section,  capacitor  input  model  with  a 
shorted  termination,   the  equation  for  which  the  roots  are  sought  is 
(from  equations  B-39): 

sinh(m+l )X  -  sinh  mX  =  0 
sinh  X[B(m+l)  -  B(m)]  =  0. 
The  case  X  =  0  has  no  meaning  for  this  model   and  so 

D(m)  =  B(m+1)  -  B(m)  =  0. 

Following  a  procedure  identical  to  that  in  Appendix  C.2.1,  the 
coefficients  of  the  "D"  polynomial  are  found  to  be 

dm(m)  =  bm(m)  (C-25) 

d.(m)  =  [b.(m)-b.(m-l)].  i=l ,2,. . . .m-1     (C-26) 

Note  that  equations  (C-25)  and  (C-26)  appear  identical  to  equations 
(C-24)  and  (C-23)  respectively,  but  that  in  the  case  of  equations 
(C-25)  and  (C-26)  the  coefficients  pertain  to  a  higher  power  of  the 
variable  C.  Equations  (C-15)  through  (C-20)  and  equations  (C-25) 
and  (C-26)  were  implemented  in  Fortran  and  solved  with  the  computer 
as  shown  in  Appendix  0.2. 


148 


C.2.3     Polynomial   for  L-Section  Model,  Open  Termination 

For  the  case  of  the  open  circuited  L-section  model,  the  equations 
for  which  the  roots  are  sought  is    (from  equation  B-46): 

2  sinh  mA  -  sinh(m+l)A  -  sinh(m-l)A  =  0 

sinh  mA(l   -  cosh  A)  =  0 

sinh  A  B(m)[l-C]  =  0. 
Again  the  case  A  =  0  adds  no  solutions,  so  that  division  by  sinh  A  can 
be  performed  assuming  A  f   0.  Finally  then 

U(m)  =  B(m)(l-C)  =  0 

D(m)  =  [bm(m)Cm"1  +  b^  (m)Cm"2  +....b2(m)  +  b^mJJO-C). 

Collecting  coefficients  of  like  powers  of  C: 

D(m)  =  -  bm(m)Cm  +  [bm(m)  -  bm  -.(mjjc"1'1  + 

m  *-  m     '         m-lv/J 

+  [b2(m)  -  b^mJjC  +  b^m)   . 

The  coefficients,  d. ,  are  then  seen  to  be 

^(m)  =  b^m)  (C-27) 

d^m)  =   [b.(m)  -  b^On)]       i=2,3,....m  (C-28) 

dm+1(m)  =  -  bm(m)  when  i=m+l  (C-29) 

Equations   (C-15)  through   (C-20)  and   (C-27)  through   (C-29)  were 
implemented  in  Fortran  and  solved  with  the  computer  as  shown  in 
Appendix  U.3. 
C.2.3     Polynomial    for  n-Section  Model,  Shorted 

For  the  case  of  the  shorted  Tt-section  model,  the  equation  for 
which  roots  are  sought  is    (from  equation  B-64): 

cosh  mA  sinh  A  =  0. 


149 


The  case  A  =  0  has  no  meaning  for  this  model  so  that 

D(m)  =  cosh  mA, 
D(m)  =  A(m), 
d. (m)  =  ai (m)    i  =  1 ,2, m+1 

This  case  is  particularly  simple  and  is  shown  implemented  in  Appendix  D.4, 
C.2.4  Polynomial  for  TT-Section  Model,  Open 

For  the  case  of  the  open  7T-section  model,  the  equation  for  which 
roots  are  sought  is  (from  equation  B-70): 

sinh  A  sinh  mA  =  0 
sinh2A  B(m)  =  0 
(1-cosh2)  B(m)  =  0 

0(m)  =  (1-C2)  B(m)  =  0  .  (C-30) 

Expanding  equation  (C-30)  and  collecting  like  powers  of  C: 

D(m)  =  bm(m)Cm+1  +  B^^mjc"1  +  Cbm.2(m)  -  bm(m)]Cm-1  +... 

...,[b2(m)  -  b4(m)]C3  +  [b^m)  -  b3(m)]C2  -  b2(m)C  -  b^m). 

The  coefficients  d.   are  then  identified  as 

d.(m)  =  -  b^m)       if     i   =1,2  (C-31) 

d.(m)  =  [b._2(m)-b.(m)]         i=3,4,...m  (C-32) 

d.(m)   =  b     Am)  i=m+l  (C-33) 

i  m-1 

d.(m)  =  b   (m)  i=m+2  (C-34) 

Equations  (C-15)  through  (C-20)  and  (C-31)  through  (C-34)  were 

implemented  in  Fortran  and  solved  with  the  computer  as  shown  in 

Appendix  0.5. 
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